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Abstract 

\  , 

The  objective  of  this  thesis  is  to  develop  a  digital  controller  using  Quanti¬ 
tative  Feedback  Theory  for  a  fighter  aircraft  with  unstable,  nonminimum  phase 
dynamics  that  meets  performance  specifications  despite  surface  failures.  Aircraft 
design  trends  for  highly  manueverable  fighter  aircraft  are  relaxing  stability  require¬ 
ments  in  order  to  increase  performance  in  the  transonic  and  supersonic  regions. 
However,  as  a  result,  the  aircraft  is  statically  unstable  in  the  subsonic  region  which 
makes  the  flight  control  system  critical  to  flight  safety.  The  conventional  approach 
to  the  flight  safety  problem  is  to  provide  multiple  redundancy  throughout  the  flight 
control  system.  However,  QFT  provides  an  alternative  to  excessive  hardware. 


The  three  controlled  states  are  the  pitch,  roll,  and  yaw  rates.  A  weighting 

matrix  is  derived  which  linearly  combines  the  nine  control  surfaces  into  three  con¬ 
i',  rlt-'vr 

trol  inpfits.  The  plant  is  converted  to  the-w>  plane  using  the  Hoffman  algorithm. 
Three  constant  gain  controllers  and  three  prefilters  are  designed  for  a  single  flight 
condition  of  0.9  Mach  and  30000  ft  altitude.  The  controllers  and  prefilters  are 
transformed  to  the  z  plane  for  simulation  purposes.  The  design  is  simulated  with 
a  healthy  plant  and  24  combinations  of  surface  failures.-The  failure  cases  consisted 
of  single  and  double  surface  impairments.  The  nonlinear  effects  of  rate  and  de¬ 
flection  limiters  on  the  actuator  models  are  added  to  the  simulations.  The  first 
set  of  simulations  consist  of  single  channel  inputs  to  determine  the  level  of  inertial 
decoupling.  A  coordinated  turn  with  55  degrees  of  bank  is  simulated  next.  The 
first  order  actuators  are  replaced  with  fourth  order  actuators  and  a  digital  filter 
is  cascaded  with  the  digital  controllers.  The  loop  transmission  is  nearly  restored 
and  a  2  degree/second  pitch  rate  command  and  90  degree/second  roll  command 
are  simulated.  An  attempt  to  extend  the  design  to  a  second  flight  condition  of  0.9 
Mach  and  200  ft  is  unsuccessful. 
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The  design  is  shown  to  meet  design  objectives  in  spite  of  surface  failures 
for  all  failure  conditions  in  the  longitudinal  mode.  Two  failure  conditions  axe  not 
satisfactorily  controlled  during  the  roll  rate  commands.  The  replacement  of  the 
first  order  actuators  with  the  fourth  order  actuators  is  shown  to  be  most  successful 
in  the  longitudinal  mode.  Excessive  oscillatory  behavior  is  found  in  the  lateral 
modes.  Quantitative  Feedback  Theory  is  shown  to  be  applicable  for  sampled-data, 
multi-input,  multi-output,  unstable,  and  nonminimum  phase  aircraft. 
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A  DIGITAL  RATE  CONTROLLER  FOR  THE  CONTROL 
RECONFIGURABLE  COMBAT  AIRCRAFT  DESIGNED  USING 
QUANTITATIVE  FEEDBACK  THEORY 


1.  Flight  Control  System  Design 

Introduction 

The  aircraft  flight  control  system  (FCS)  has  evolved  into  a  highly  redun¬ 
dant,  flight-critical,  and  maintenance  intensive  system  at  the  expense  of  aircraft 
performance  5  .  On  current  high  performance  aircraft,  the  conventional  FCS  de¬ 
sign  approach  has  resulted  in  large,  high  authority  control  aerosurfaces  to  generate 
the  rotational  moments  needed  for  both  performance  and  flight  safety  [5\  Failure 
of  any  of  these  surfaces  can  result  in  a  loss  of  control  in  the  affected  dimension 
and  probable  loss  of  the  aircraft.  Air  Force  operational  policy  is  a  result  of  the 
criticality  of  the  FCS;  no  aircraft  can  start  a  mission  without  a  100%  operational 
FCS  and  the  FCS  must  be  able  to  withstand  two  failures  during  a  mission  without 
any  loss  of  capability  5  .  The  costs  imposed  by  such  requirements  on  the  FCS 
are  weight,  complexity,  and  reliability  which  directly  affects  payload,  range,  and 
performance  [5’.  If  the  value  of  any  one  aerosurface  in  terms  of  flight  safety  could 
be  reduced  without  a  loss  in  performance,  then  the  damage  of  that  surface  might 
be  sustained  without  the  potential  loss  of  the  mission  or  the  aircraft.  One  method 
of  reducing  the  value  of  a  control  surface  is  to  design  an  FCS  which  is  tolerant  of 
control  surface  damage.  In  other  words,  the  FCS  design  should  continue  to  meet 
the  performance  specifications  despite  an  uncertainty  in  the  available  control  force 
from  a  control  surface.  Such  a  reduction  in  the  criticality  of  a  control  surface  would 


allow  the  simplification  of  the  design  of  many  of  the  components  of  the  FCS  with 
a  parallel  increase  in  mission  effectiveness. 

Quantitative  Feedback  Theory  (QFT)  is  a  design  approach  which  specifically 
confronts  the  problem  of  uncertainty  in  feedback  systems  in  a  manner  which  en¬ 
courages  the  explicit  trade-off  of  design  parameters.  The  philosophy  behind  QFT 
is  firmly  rooted  in  the  classical  frequency  approach  to  feedback  design.  As  Pro¬ 
fessor  Horowitz  points  out  in  his  book,  Synthesis  of  Feedback  Systems ,  the  use  of 
feedback  is  motivated  by  uncertainty  of  the  plant  dynamics,  plant  parameters,  and 
plant  inputs  in  the  form  of  disturbances.  This  uncertainty  causes  the  desired  plant 
output  and  actual  plant  output  to  be  different.  Feedback  is  employed  to  compare 
the  intended  output  with  the  actual  output  and  adjust  the  input  to  the  plant  ac¬ 
cordingly  1 10,.  The  difficulty  lies  in  determining  how  much  feedback  to  employ; 
too  little  feedback  and  the  uncertainty  is  not  adequately  controlled  and  too  much 
feedback  and  the  loop  bandwidth  increases  to  an  unacceptable  level.  However, 
QFT  optimizes  the  amount  of  feedback  for  a  given  amount  of  parameter  uncer¬ 
tainty  [11,12,13,14,16:.  The  optimization  process  is  accomplished  by  examining 
the  costs  and  benefits  of  feedback  across  the  frequency  spectrum  of  interest. 

One  of  the  first  decisions  which  must  be  made  in  the  design  of  a  control 
system  is  to  determine  the  amount  of  uncertainty  which  will  be  controlled  [17,19^ 
The  level  of  feedback  required  increases  with  increasing  uncertainty.  Eventually, 
the  amount  of  controllable  uncertainty  is  limited  by  the  maximum  permissible 
loop  bandwidth.  This  example  is  typical  of  a  control  system  design;  the  continual 
balancing  of  competing  requirements.  The  advantage  of  the  QFT  approach  is  that 
many  of  the  trade-offs  are  explicit.  The  design  choices  are  made  with  the  impact  on 
design  parameters  such  as  bandwidth,  noise  susceptibility,  or  phase  margin  clearly 
defined.  This  transparency  makes  QFT  ideally  suited  for  designing  an  aircraft 
FCS. 
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Problem  Description 

The  Control  Reconfigurable  Combat  Aircraft  (CRCA)  project  at  the  Air 
Force  Wright  Aeronautical  Laboratory  (AFWAL)  Flight  Dynamics  Laboratory 
(FDL)  is  an  ongoing  effort  to  integrate  a  robust  flight  control  system  with  a  failure 
detection  system.  The  aircraft  model  is  a  hypothetical  advanced  tactical  aircraft 
modeled  in  a  six  degree-of-freedom  nonlinear  simulation.  The  existing  flight  con¬ 
trol  system  is  a  proportional  plus  integral  (PI)  controller  in  the  feed-forward  loop 
for  both  the  longitudinal  and  lateral  axes  ;26i.  While  the  current  FCS  meets  or 
exceeds  most  of  the  requirements  for  a  category  II  (fighter)  aircraft,  several  areas 
for  improvement  exist.  Specifically,  unnecessary  pitch  rate  overshoot  is  evident 
and  yaw  and  roll  rates  are  coupled  under  certain  failure  conditions.  The  purpose 
of  this  thesis  is  to  replace  the  existing  flight  control  system  with  one  designed  using 
QFT  without  a  reliance  on  a  failure  detection  system.  The  primary  goal  for  this 
thesis  is  to  minimize  pitch  rate  overshoot  and  decouple  yaw  and  roll  rates  in  spite 
of  control  surface  failure.  A  secondary  goal  is  to  demonstrate  the  applicability  of 
QFT  to  an  unstable,  nonminimum  phase,  sampled-data  flight  control  system. 

Background 

Advanced  fighter  design  concepts  are  relaxing  the  static  stability  require¬ 
ments  in  both  the  longitudinal  and  lateral  axes  to  improve  transonic  and  super¬ 
sonic  performance  [4,26’.  The  F-16,  the  AFTI  F-16,  and  the  X-29  are  examples  of 
aircraft  which  are  unstable  in  the  longitudinal  axis.  In  addition,  dynamic  analysis 
of  modern  fighter  aircraft  reveals  that  unacceptable  levels  of  inertial  cross-coupling 
can  occur  during  certain  manuevers  such  as  high-speed  rolls  [7j.  The  combination 
can  place  such  large  demands  on  the  control  system  that  the  FCS  becomes  the 
limiting  factor  on  aircraft  performance.  For  instance,  enhancing  transonic  perfor¬ 
mance  by  moving  the  center  of  gravity  backwards  is  limited  by  the  necessity  of 
maintaining  sufficient  nose-down  pitching  moment  in  high  angle  of  attack  manuev- 
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ers  to  prevent  departure.  In  addition  to  maintaining  control  during  a  commanded 
input,  sufficient  control  power  must  also  exist  at  all  points  in  the  flight  envelope 
to  counter  the  effects  of  inertial  cross-coupling,  airframe  instabilities,  and  turbu¬ 
lence  [4!.  The  control  surfaces  can  be  increased  in  size  to  generate  larger  rotational 
moments,  but  with  the  attendant  cost  of  increasing  the  criticality  of  the  aerosur- 
faces.  If  the  surface  size  is  increased  to  provide  the  control  power  necessary  to 
maintain  stability,  then  the  loss  of  that  surface  will  precipitate  the  loss  of  control 
of  the  aircraft. 

The  objective  of  a  FCS  design  for  a  statically  unstable  aircraft  is  to  provide 
adequate  control  power  to  meet  the  performance  requirements  while  minimizing 
overdesign  of  the  FCS  without  compromising  flight  safety.  The  increased  perfor¬ 
mance  capabilities  must  be  balanced  against  the  adverse  effects  caused  by  stabi¬ 
lizing  the  divergent  airframe  roots  [4  .  Obviously,  a  10  g  aircraft  which  cannot 
adequately  reject  disturbance  inputs  is  not  a  high  performance  aircraft.  In  ad¬ 
dition,  the  design  must  also  take  into  account  the  performance  limitations  of  the 
hardware  and  uncertainty  in  the  plant  parameters.  A  design  which  requires  un¬ 
realistic  control  surface  deflections  to  meet  performance  requirements  or  does  not 
control  the  plant  throughout  the  flight  envelope  is  obviously  unacceptable.  Perfor¬ 
mance  vs  disturbance  rejection  and  performance  vs  hardware  limitations  are  two 
examples  of  the  kinds  of  trade-offs  faced  in  the  design  of  a  FCS.  The  difficulty 
lies  in  relating  these  trade-offs  to  design  parameters  such  as  bandwidth  and  how  a 
particular  design  step  impacts  the  outcome  of  a  trade-off. 

A  flight  control  system  must  manage  numerous  dynamic  uncertainties.  The 
healthy  aircraft  is  subjected  to  plant  parameter  variation  due  to  dynamic  loading 
and  noise  inputs  due  to  gusts.  In  addition,  partially  understood  structural  modes, 
sensor  noises,  and  unknown  effects  due  to  aging  further  complicate  the  design  pro¬ 
cess.  Finally,  an  aircraft  with  damaged  control  surfaces  is  subjected  to  even  further 
plant  parameter  variations.  One  approach  to  the  design  of  a  FCS  is  to  quantify  the 
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amount  of  uncertainty  before  the  design  process  begins.  For  instance,  the  dynamic 
effects  of  a  missing  surface  can  be  predicted  with  wind  tunnel  data  derived  from 
a  model  with  the  surface  removed.  Once  a  set  of  dynamic  equations  describing 
a  range  of  plant  variations  is  derived,  a  controller  which  is  tolerant  to  the  set  of 
modeled  impairments  can  be  designed.  At  the  same  time,  equal  importance  should 
be  placed  on  reducing  the  amount  of  overdesign  in  the  controller.  Overdesign  is 
defined  as  providing  more  feedback  than  necessary  to  reduce  the  plant  sensitivity 
to  a  parameter  variation  1 10  .  If  the  flight  control  system  is  robust  to  parameter 
variations  and  has  a  minimal  amount  of  overdesign,  the  overall  complexity  of  the 
flight  control  system  could  be  reduced.  The  objectives  of  the  CRCA  program  are  to 
increase  reliability,  ease  maintainability,  and  increase  weapon  effectiveness.  These 
objectives  are  attainable  at  least  in  part  by  using  a  flight  control  system  designed 
using  QFT. 

The  loop  bandwidth  of  flight  control  system  is  one  measure  of  how  well  the 
desired  objectives  are  balanced  against  the  imposed  costs  of  feedback.  The  band¬ 
width  defines  the  frequency  range  where  the  open-loop  gain  must  be  greater  than 
one  in  order  to  obtain  good  system  performance  while  maintaining  tractable  dis¬ 
turbance  rejection  characteristics  j4].  As  the  bandwidth  is  increased,  the  amount 
of  disturbance  rejection  increases.  However, the  response  of  the  plant  to  high  fre¬ 
quency  components  such  as  sensor  noise  and  structural  modes  is  also  increasing. 
If  the  bandwidth  becomes  too  large,  the  control  system  can  introduce  high  fre¬ 
quency  flutter  to  the  aerosurfaces  which  can  lead  to  saturation  of  the  actuators. 
The  obvious  trade-off  is  to  maximize  bandwidth  until  the  frequency  range  of  the 
sensor  noise  or  structural  modes  are  encountered.  The  obstacle  for  many  design 
approaches  lies  in  finding  precise  relationships  between  the  desired  design  mea¬ 
sures  of  merit  and  specific  design  steps.  QFT  on  the  other  hand  is  particularly 
well  suited  to  explicitly  describing  the  relationship  between  frequency  design  spec¬ 
ifications  such  as  gain  margin  or  phase  margin  and  controller  design.  The  effect 
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of  adding  a  pole  or  zero  to  the  controller  on  bandwidth,  phase  or  gain  margin,  or 
crossover  frequency  is  graphically  portrayed  when  the  loop  transmission  is  plotted 
on  the  Nichols  chart.  Trade-offs  can  be  made  with  the  benefits  and  costs  explicitly 
shown. 

One  example  of  how  QFT  is  ideally  suited  to  FCS  design  is  demonstrated 
by  displaying  the  limits  on  the  amount  of  uncertainty  that  the  use  of  feedback 
can  control.  Templates  are  used  in  QFT  to  graphically  portray  the  phase  angle 
and  magnitude  characteristics  of  a  transfer  function  with  a  range  of  uncertainty  at 
a  particular  frequency.  The  larger  the  uncertainty,  the  larger  the  template  tends 
to  become.  The  templates  are  used  to  determine  the  stability  and  performance 
bounds  which  guide  the  shaping  of  the  loop  transmission  function.  As  the  size  of 
the  templates  increase,  the  resultant  bounds  force  the  frequency  of  the  crossover 
point  to  increase.  The  open  loop  crossover  point  is  a  rough  approximation  of 
where  the  dynamics  of  the  closed  loop  system  starts  to  attenuate  rather  than 
amplify  information  from  the  higher  frequencies.  If  the  bandwidth  extends  into 
the  regions  where  noise  or  structural  modes  dominate,  then  the  costs  of  feedback 
have  overwhelmed  the  desired  benefits.  Therefore,  the  templates  allow  the  designer 
to  set  realistic  performance  goals  graphically  revealing  the  limitations  of  feedback. 

A  second  example  of  how'  QFT  is  applicable  to  the  design  of  flight  control 
systems  is  in  the  area  of  stability  margins.  The  military  design  specification  mil 
std  9490D  requires  a  high  frequency  gain  margin  of  6  dB.  The  stability  margins  are 
built  into  the  design  to  cope  with  unmodeled  aspects  of  the  plant  such  as  higher 
order  dynamics  or  sensor  noise  [4j.  If  the  design  technique  employed  does  not  allow 
visibility  into  characteristics  such  as  the  loop  transmission  frequency  response,  the 
compensator  may  meet  performance  specifications  but  might  be  overdesigned.  In 
such  a  case,  the  frequency  response  of  the  loop  transmission  function  will  needlessly 
extend  into  the  higher  frequencies.  The  result  can  be  excessive  high  frequency  noise 
amplification  and  potential  plant  saturation.  QFT  avoids  this  situation  during  the 
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loop  shaping  process.  Loop  shaping  is  performed  by  plotting  the  loop  transmis¬ 
sion  function  on  a  Nichols  chart.  The  performance  specifications  are  embedded  in 
design  equations  which  are  used  to  plot  boundaries  on  the  Nichols  chart.  The  com¬ 
pensator  is  designed  to  meet  the  performance  boundaries  within  the  constraints  of 
the  stability  margins.  A  choice  for  the  compensator  is  quickly  shown  to  be  ade¬ 
quate  based  on  the  loop  transmission  plot  relative  to  the  performance  bounds.  As 
the  compensator  is  developed,  critical  design  parameters  such  as  gain  margin  or 
susceptibility  to  noise  amplification  are  continually  displayed  on  the  chart.  Overde¬ 
sign  is  avoided  since  the  amount  of  feedback  can  be  fine-tuned  with  respect  to  the 
design  parameters.  The  “guess  and  simulater  approach  is  unnecessary  and  design 
time  is  greatly  reduced. 

I 

Approach 

The  development  of  this  thesis  is  broken  into  several  parts.  First,  the  basic 
*  plant  models  for  the  nominal  case  and  the  failure  set  are  derived  from  aerodynamic 

derivatives.  Second,  a  weighting  matrix  is  derived  which  linearly  combines  the  nine 
columns  of  the  input  matrix  into  three  columns.  Third,  the  continuous  plant  and 
actuators  are  transformed  to  the  discrete  wt  plane.  Fourth,  the  three  compensators 
r  are  designed  in  the  wt  plane  using  the  Quantitative  Feedback  Theory  approach. 

Finally,  the  closed  loop  system  is  simulated  in  the  z  plane  using  the  dynamic  system 
modeling  capabilities  of  MatrixX. 

Assumptions 

The  following  assumptions  form  the  basis  from  which  the  work  performed  in 
i  this  thesis  is  begun. 

1.  The  aerodynamic  data  supplied  by  AFWAL/FDL  adequately  represents  the 


aircraft. 


2.  The  software  provided  by  AFW’AL  FDL  used  to  obtain  the  linear  aircraft 
state  space  model  operates  correctly. 

3.  The  fixed  sampling  rate  of  80  Hz  is  sufficient. 

4.  The  mass  of  the  aircraft  remains  constant  throughout  the  simulation  time 
period. 

5.  The  aircraft  engine  thrust  remains  constant  through  the  simulation  time 
period. 

6.  The  aircraft  can  be  treated  as  a  rigid  body  for  the  purposes  of  deriving  the 
linear  equations  of  motion. 

7.  The  earth  surface  is  the  inertial  reference  frame. 

8.  The  atmosphere  is  fixed  to  the  surface  of  the  earth. 

9.  No  failure  detection  schemes  will  be  employed  to  augment  the  constant  gain 
controller  to  maintain  stability. 

10.  The  actuators  have  no  uncertainty  across  the  frequency  spectrum. 

11.  Effects  of  finite  wordlength  will  not  be  investigated. 

Presentation 

This  thesis  is  divided  into  six  chapters.  The  development  of  the  basic  aircraft 
model  in  the  s  plane  and  the  weighting  matrix  derivation  are  presented  in  chapter 
2.  The  third  chapter  presents  the  effective  plant  development  in  the  wt  plane  and 
the  digital  QFT  design  procedures.  Chapter  4  presents  the  design  of  the  digital 
compensators  and  prefilters.  Simulation  results  are  presented  in  chapter  5.  Finally, 
conclusions  and  recommendations  are  presented  in  the  sixth  chapter. 


II.  Model  Development 


Overview 

The  objective  of  this  chapter  is  to  outline  the  procedures  used  to  develop 
the  basic  plant  transfer  function  matrix,  the  weighting  matrix,  and  the  effective 
plant  transfer  function  matrix  in  the  continuous  s  plane.  The  basic  plant  is  defined 
as  the  8x9  transfer  function  matrix.  The  effective  plant  is  defined  as  basic  plant 
postmultiplied  by  the  weighting  matrix.  Three  states  are  used  in  feedback  so  the 
final  transfer  function  matrix  is  a  square  3x3  matrix.  A  linear  state  space  model 
of  the  aircraft  is  derived.  The  equivalent  basic  plant  transfer  function  matrix  is 
obtained  in  the  continuous  s  plane.  A  weighting  matrix  is  developed  which  linearly 
combines  the  nine  columns  of  the  input  matrix  into  three  columns.  The  effective 
plant  is  derived  in  the  continuous  plane  by  postmultiplication  of  the  plant  matrix 
with  the  weighting  matrix.  The  actuators  are  cascaded  with  the  effective  plant 
model.  The  effective  plant  is  transformed  to  the  wi  plane  and  the  determinant 
of  the  3x3  transfer  function  matrix  is  taken.  The  effective  plant  is  shown  to  be 
minimum  phase  for  25  of  the  potential  28  failure  cases  for  one  aerodynamic  flight 
condition  and  for  16  of  17  potential  failure  cases  for  a  second  aerodynamic  flight 
condition. 

Aircraft  Description 

The  CRCA  is  a  modified  version  of  a  NASA/Grumman  advanced  tactical 
fighter  (ATF)  class  aircraft  model.  The  CRCA  is  a  “paper  airplane”;  a  math 
model  simulated  on  the  Air  Force  Simulation/Rapid-Prototyping  Facility  located 
at  the  Air  Force  Wright  Aeronautical  Laboratories  (AFWAL)  Flight  Dynamics 
Laboratory  (FDL).  The  aircraft  is  representative  of  a  1995  fighter  type  aircraft. 
At  low  airspeeds,  the  aircraft  is  12%  unstable  in  the  longitudinal  mode.  There  are 
nine  control  surfaces:  two  canards  with  thirty  degrees  of  dihedral,  four  trailing-edge 
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flaps,  two  elevators,  and  a  rudder.  The  canards  provide  the  main  pitch  moment 
and  a  smaller  amount  of  yaw  and  roll  moments.  The  two  trailing-edge  flaps  on 
each  wing  are  used  to  provide  both  rolling  and  yawing  moments.  The  two  elevators 
provide  both  pitch  and  roll  moments.  The  rudder  supplies  yaw  moment  26  . 


Figure  2.1.  Control  Reconfigurable  Combat  Aircraft 


The  design  for  this  thesis  is  developed  using  data  provided  for  the  Air  Combat 
■  Manuever  (ACM)  Entry  condition.  The  condition  represents  the  aerodynamic 

effects  of  a  fighter  maneuvering  at  0.9  Mach  at  30000  ft.  Combat  maneuvering 
requires  the  aircraft  to  roll  rapidly  without  inducing  a  loss  of  control  due  to  inertial 
cross-coupling,  to  quickly  induce  maximum  g- loading,  and  to  maintain  coordination 
to  reduce  drag.  The  initial  aircraft  orientation  is  the  trim  condition  with  the  wings 
level.  After  completion  of  the  design,  the  unmodified  controller  is  also  simulated  in 
the  Terrain-Following,  Terrain-Avoidance  (TFTA)  condition.  The  two  conditions 
are  similar  in  terms  of  g  loading  but  the  Q  factor  (dynamic  loading)  for  the  TFTA 
manuever  is  much  larger.  The  TFTA  represents  the  effects  of  flying  at  200  ft  at 
0.9  Mach  at  low  angles  of  attack  with  rapid  but  limited  changes  in  pitch  angle  and 
)  bank  angle  [26'. 
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Basic  Plant 


The  trim  or  nominal  case  and  subsequent  failure  cases  are  developed  from 
nondimensional  aerodynamic  derivatives.  The  derivatives  axe  the  result  of  a  non¬ 
linear  simulation  of  the  aircraft  drawing  on  a  database  made  up  of  a  blend  of 
wind  tunnel  data  and  the  estimated  aerodynamic  effects  of  partial  aerosurface 
losses.  The  aircraft  surfaces  were  trimmed  to  balance  the  resulting  moments  and 
the  thrust  was  adjusted  to  compensate  for  any  change  in  drag  to  maintain  the 
initial  velocity.  The  nonlinear  equations  of  motion  were  then  perturbed  about  an 
equilibrium  point  to  generate  the  nondimensional  aerodynamic  derivatives.  The 
derivatives  were  supplied  to  the  author  by  AFWAL/FDCL. 

The  linear  equations  of  motion  are  developed  using  two  programs  provided 
by  AFWAL/FDCL  {3}  The  first  program  uses  the  nondimensional  coefficients  as 
input  data  to  derive  the  dimensional  coefficients.  The  second  program  uses  the 
dimensional  coefficients  to  develop  the  linear  state  equations.  The  linear  equations 
of  motion  are  developed  in  the  standard  state  space  format  of 

x  =  Ax  4-  Bu 

y  —  Cx  +  Du  (2.1) 

where  x  is  an  9x1  state  variable  vector, 
u  is  a  9x1  state  input  vector. 

A  is  an  9x9  state  dynamics  matrix. 

B  is  an  9x9  state  input  matrix. 

C  is  a  3x9  state  output  matrix. 

D  is  a  3x9  state  feed  forward  matrix. 
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The  state  variable  vector  for  the  CRCA  is  given  by 


x  = 


u 

w 

Q 

e 

<f> 

0 

p 

r 


The  state  control  input  vector  for  the  CRCA  is  given  by 


u  = 


left  canard 
right  canard 
left  outer  flap 
left  inner  flap 
right  outer  flap 
right  inner  flap 
left  elevator 
right  elevator 
rudder 


(2.2) 


(2.3) 


The  linear  state  space  model  consists  of  a  nine  state  A  matrix  and  a  nine 
input  B  matrix.  The  A  matrix  is  block  diagonal  indicating  the  longitudinal  and 
lateral  modes  are  decoupled.  The  B  matrix  is  not  diagonal;  most  control  inputs 
affect  both  the  longitudinal  and  lateral  states;  the  only  single  axis  surface  is  the 
rudder.  The  A  matrix  is  rank  deficient  so  the  equation 


xfj  =  r  (2.4) 

is  removed  by  deleting  the  fifth  row  and  fifth  column  of  the  A  matrix  and  the  fifth 
row  of  the  B  matrix.  The  A  matrix  is  now  full  rank  and  therefore  invertible.  The 
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A  matrix  is  now  8x8,  the  B  matrix  is  8x9.  and  the  state  vector  is  8x1.  The  A 
matrix  for  the  nominal  condition  of  the  ACM  manuever  is  given  by: 
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The  B  matrix  for  the  nominal  condition  is  given  by: 
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-0.111 
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-0.049 

0.003 

0.002 

-0.003 

-.002 

0.002 
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(2.6) 


The  three  state  variables  which  are  controlled  by  feedback  are  roll  rate,  pitch 
rate,  and  yaw  rate.  The  required  C  matrix  is  then  given  by: 


0  0  1.0  0  0  0  0  0 
0  0  0  0  0  0  1.0  0 
0  0  0  0  0  0  0  1.0 


(2.7) 


Since  there  are  no  feed-forward  terms  the  D  matrix  is  a  3x9  zero  matrix. 


Using  the  expression, 

=  <?(«)<*/  -  +  D(.) 

V[s) 


(2.8) 
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Aerosurface  Deflection  Limit  Rate  Limit 


Canards  —60  degrees  -30  degrees  100  degrees /second 


Trailing  Edge  Flaps  ;  -30  degrees  ;  -30  degrees  100  degrees /second 
Rudder  j  -t-20  degrees  i  -20  degrees  100  degrees /second 


Table  2.1.  Deflection  and  Rate  Limits 

the  corresponding  transfer  functions  for  the  basic  nine-input,  three-output  plant 
are  derived  6  .  The  basic  plant  transfer  functions  are  in  appendix  C.  The  transfer 
functions  for  the  commanded  pitch  rate  have  an  unstable  root  at  s  =  1.694.  The 
commanded  roll  rate  transfer  functions  are  stable  with  a  small  nonminimum  phase 
zero.  The  commanded  yaw  rate  transfer  functions  are  also  stable  but  with  two 
right  half  plane  complex  zeros. 

The  actuators  for  the  CRCA  are  modeled  using  two  fourth  order  approxima¬ 
tions.  The  actuators  for  the  canards  and  flaps  are  modeled  by: 

Y(s)  =  _ 1-1816 -107 _  .  > 

U{s)  (s  ^  20.20)(s  +  52.55  ±j48.34)(s+  114.7)  1  '  ’ 

The  actuator  for  the  rudder  is  modeled  by: 


Y(«)  _  1.1551  107 

U(s)  ~  (s  -  20.6)(s  -  54.1  =  j46.l)(s  +  111) 


(2.10) 


The  deflection  and  rate  limits  for  the  control  surfaces  are  supplied  by  AFWAL / FDCL 
and  are  given  in  Table  2.1. 

A  linear  simulation  of  the  step  response  of  both  models  is  performed.  The 
results  in  Figure  2.2  show  that  the  two  models  have  virtually  identical  time  re¬ 
sponses  to  step  inputs.  Therefore,  the  rudder  actuator  model  is  replaced  by  the 
canard  or  flaperon  actuator  model  with  no  loss  of  fidelity.  The  actuator  model 
is  cascaded  with  each  basic  plant  transfer  function  and  is  considered  to  have  no 
uncertainty  throughout  the  frequency  range. 
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Figure  2.2.  Time  Response  of  Canard /Flaperon  and  Rudder  Actuator  Models 


Failure  Cases 

The  aerodynamic  data  for  27  simulated  single  and  double  surface  failures 
for  the  ACM  entry  manuever  and  17  simulated  failures  for  the  TFTA  manuever 
was  supplied  by  the  AFWAL/FDCL.  The  data  is  used  to  derive  the  corresponding 
state  space  models  using  the  same  approach  described  for  the  nominal  case.  The 
complete  set  for  each  aerodynamic  condition  consisting  of  the  trimmed,  unimpaired 
case  and  the  simulated  surface  failures  is  referred  to  in  this  thesis  as  the  design 
set.  The  simulated  failures  are  loss-of-effectiveness  failures  and  fail-to-trail  failures 
of  single  and  double  surfaces.  The  loss-of-effectiveness  failures  simulate  battle 
damage  which  removes  25%,  50%,  and  100%  of  the  affected  aerosurface.  The  fail- 
to-trail  failures  simulate  a  surface  locked  at  0.7  times  the  local  angle  of  attack.  The 
complete  list  is  located  in  appendix  B. 

The  transfer  functions  derived  from  the  design  set  for  ACM  entry  are  used  to 
develop  the  templates.  The  templates  are  used  to  develop  the  stability  boundaries. 
Once  a  control  system  is  designed  which  meets  the  stability  boundaries,  the  closed 


loop  system  is  simulated  with  the  complete  design  set.  The  objective  is  to  examine 
the  ability  of  the  constant  gain  QFT  controller  to  continue  meeting  performance 
goals  in  spite  of  parameter  variation  due  to  surface  failures. 

Development  of  the  Weighting  Matrix 

A  weighting  matrix  is  developed  which  converts  the  basic  plant  transfer  func¬ 
tion  matrix  into  a  3x3  effective  plant  transfer  function  matrix.  The  requirement  for 
a  square  transfer  function  matrix  is  due  to  the  manner  in  which  QFT  approaches 
the  multi-input,  multi-output  (MIMO)  design  problem.  An  integral  step  of  the  de¬ 
sign  procedure  is  to  convert  a  nxn  plant  into  n2  multi-input,  single-output  (MISO) 
plants.  The  conversion  process  requires  the  inversion  of  the  MIMO  plant  matrix. 
Therefore,  the  MIMO  transfer  function  matrix  must  be  square.  The  weighting 
matrix  must  accomplish  two  unrelated  objectives.  First,  the  weighting  matrix 
combines  the  terms  of  the  control  input  matrix  so  that  the  determinant  of  the  ef¬ 
fective  plant  matrix  is  minimum  phase  for  each  member  of  the  design  set.  Second, 
the  weighting  matrix  spreads  the  control  authority  required  to  achieve  the  desired 
output  across  multiple  control  surfaces.  Both  goals  are  accomplished  by  linearly 
combining  the  nine  columns  of  the  input  control  matrix  into  three  columns  and 
thereby  creating  a  three  input  control  matrix. 

The  first  objective  for  the  weighting  matrix  is  necessary  in  order  to  meet  the 
performance  specifications.  The  beneficial  effects  of  feedback  such  as  disturbance 
rejection  are  greatly  reduced  if  the  MIMO  plant  is  not  minimum  phase.  An  MIMO 
plant  is  defined  to  be  minimum  phase  if  the  determinant  of  the  plant  numerator 
matrix  has  no  roots  in  the  right-half  plane  [10].  The  explanation  for  the  reduc¬ 
tion  in  performance  of  a  nonminimum  phase  plant  lies  in  the  phase  response  of  a 
transfer  function  due  the  nonminimum  phase  zero.  At  arbitrarily  low  frequencies,  a 
minimum  phase  simple  zero  makes  a  phase  contribution  of  zero  degrees,  but  a  non¬ 
minimum  phase  zero  makes  a  phase  contribution  of  180  degrees.  As  the  frequency 
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Figure  2.3.  Effective  Plant 


increases,  the  minimum  phase  zero  begins  to  increase  in  phase  until  a  maximum 
of  90  degrees  is  reached.  On  the  other  hand,  the  nonminimum  phase  zero  will 
decrease  in  phase  until  90  degrees  is  also  reached.  Therefore,  at  low  frequencies, 
the  phase  margins  required  for  stability  are  eroded  by  a  nonminimum  phase  zero 
and  consequently  the  gain  may  need  to  be  lowered  to  maintain  stability.  The  low¬ 
ered  gain  reduces  the  bandwidth  over  which  the  beneficial  effects  of  feedback  can 
operate  and  accordingly  performance  is  reduced. 

The  second  objective  for  the  weighting  matrix  is  to  disperse  the  total  required 
moment  for  a  given  input  across  multiple  control  surfaces.  The  pilot  inputs  to  the 
closed  loop  CRCA  for  this  thesis  are  pitch,  roll,  and  yaw  rate.  However,  the  actual 
control  inputs  to  the  aircraft  are  the  nine  aerosurface  movements.  The  weighting 
matrix  translates  the  desired  pilot  inputs  into  control  inputs  for  the  CRCA  by 
taking  advantage  of  the  multiaxis  response  of  most  of  the  CRCA  surfaces.  All 
surfaces  except  the  rudder  are  used  to  implement  a  pitch  or  roll  rate  command 
and  all  the  surfaces  are  used  for  the  yaw  command.  Therefore,  all  the  surfaces  are 
employed  to  implement  a  pilot  command.  The  individual  terms  of  the  weighting 


matrix  determine  the  level  of  effort  required  of  each  surface  for  each  pilot  command. 
The  magnitude  and  signs  of  the  weighting  matrix  terms  must  be  such  that  the 
resulting  aerosurface  movements  are  coordinated  and  the  deflections  are  within 
reasonable  limits  for  given  desired  output.  An  additional  feature  of  the  weighting 
matrix  is  the  criticality  of  any  particular  surface  can  be  reduced  by  spreading 
control  authority  across  multiple  surfaces.  If  a  surface  is  impaired,  the  required 
moments  can  still  be  generated  by  the  remaining  surfaces.  The  reconfigurability 
of  the  control  inputs  paths  allows  aircraft  performance  to  degrade  in  incremental 
steps  with  multiple  surface  failures. 

The  conversion  from  the  basic  plant  to  an  effective  plant  is  accomplished  by 
post-multiplying  the  basic  plant  transfer  function  matrix  by  the  weighting  matrix 
in  the  continuous  s  plane. 

=  {C(«)(»J  -  A)-'B(s)  +  Z>(«)}  •  A  (2.11) 

Since  the  D  matrix  is  a  matrix  of  zeros,  only  the  B  matrix  is  post-multiplied  by 
the  weighting  matrix. 

Techniques  for  deriving  the  terms  of  the  weighting  matrix  can  be  best  de¬ 
scribed  as  opaque.  Earlier  QFT  MIMO  thesis  efforts  relied  on  previous  experience 
and  aerodynamic  insight  to  determine  the  weightings  for  a  surface.  However,  a 
recent  thesis  developed  a  method  for  determining  the  relative  magnitudes  of  the 
weighting  matrix  elements.  The  method  is  called  Specified  Output  because  the 
weighting  matrix  is  derived  based  on  the  desired  steady  state  output  [81.  The 
expression  for  the  effective  plant,  Pe[s),  in  the  s  plane  is  given  by 

(s)  =  Pb(s)  •  A  (2.12) 

In  the  steady  state,  if  input  equals  the  output,  then  the  left  hand  side  of  the 
expression  in  equation  (2.12)  can  be  viewed  as  equivalent  to  a  3x3  identity  matrix. 
At  the  same  time,  zeros  in  the  off-diagonal  positions  indicates  the  lack  of  coupling 
between  the  commanded  state  and  the  other  two  uncommanded  states. 


If  the  product  of  the  two  terms  on  the  right  of  the  expression  in  equa¬ 
tion  (2.12)  is  an  identity  matrix,  then  delta  must  be  equal  to  the  inverse  of  Pt(s). 
But  since  Pb(s)  is  a  3x9  matrix  and  is  not  invertible,  a  pseudoinverse  must  be 
used  [8j.  The  pseudoinverse  is  defined  as 

{AM}"1  -  {a'M-AM}  '-A,'M  (213) 

Manipulating  the  Laplace  form  of  the  plant  transfer  functions  to  develop 
the  weighting  matrix  is  extremely  cumbersome.  Instead,  the  weighting  matrix  for 
this  thesis  is  developed  numerically  by  substituting  jw  for  s  in  the  plant  transfer 
function  matrices  used  in  equation  (2.13).  The  magnitude  of  the  resulting  complex 
value  of  equation  (2.13)  is  found  for  five  frequencies:  w  =  1,2,4, 8,  and  16  rad/ sec. 
The  five  magnitudes  are  summed  and  the  average  determined.  The  signs  of  the 
individual  weighting  matrix  terms  are  determined  by  examination  of  the  frequency 
response  of  the  27  transfer  functions  of  the  basic  plant  [17^.  Multiplication  of  a 
term  by  -1  is  used  to  force  like  aerosurfaces  (left  canard,  right  canard  and  so 
on)  to  have  similar  phase  characteristics  which  indicates  the  surfaces  are  working 
together.  The  sign  convention  is  overlaid  on  the  magnitude  terms.  The  result  is 
normalized  by  some  arbitrary  factor  to  avoid  heavily  weighting  the  control  surfaces 
in  terms  of  surface  deflection  rates  or  angles.  The  weighting  matrix  derived  for  this 
thesis  is 
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The  derivation  of  the  weighting  matrix  using  the  method  described  does 
not  claim  to  be  the  optimal  approach.  If  dynamics  are  used  rather  than  only  j 

magnitudes  for  the  individual  elements,  the  phase  and  magnitude  characteristics 
of  like  surfaces  could  be  matched  closely  over  a  large  frequency  range.  However,  the 
objective  of  this  thesis  is  to  derive  an  adequate  matrix  and  not  the  optimal  matrix.  j 

The  criteria  for  an  adequate  matrix  are  to  produce  a  minimum  phase  determinant 
of  the  effective  plant  for  each  member  of  the  design  set  in  the  discrete  domain. 

The  determinant  of  the  effective  plant  is  minimum  phase  if  there  are  no  right  half 

plane  zeros:  this  definition  is  slightly  modified  in  the  discrete  u’  plane  and  is  further  * 

discussed  in  the  next  chapter.  Using  the  weighting  matrix  in  equation  (2.14),  the 

determinant  of  the  effective  plant  is  shown  in  the  next  chapter  to  be  minimum 

phase  for  the  nominal  case  for  both  ACM  entry  and  TFTA  and  for  24  of  the  j 

possible  27  cases  for  the  ACM  manuever  entry  and  for  16  of  the  17  cases  for  the 

TFTA  manuever. 

j 

Performance  Specifications 

An  integral  part  of  the  QFT  approach  to  designing  a  control  system  is  to 
embed  the  performance  specifications  in  the  design  equations.  The  objective  is  to 
allow  the  designer  to  intelligently  compare  design  choices  with  the  effect  on  the 
performance  specifications  clearly  displayed.  QFT  develops  three  types  of  spec¬ 
ifications  in  the  form  of  point-wise  frequency  bounds:  stability  bounds,  tracking 
bounds,  and  disturbance  bounds.  The  stability  bounds  define  the  region  which 
can  not  be  penetrated  by  the  loop  transmission  function.  The  tracking  and  dis¬ 
turbance  bounds  are  magnitude  and  frequency  bounds  on  the  Nichols  chart  which 
the  loop  transmission  must  not  exceed  if  the  specifications  are  to  be  met.  Typi¬ 
cally,  the  tracking  bounds  dominate  in  the  low  frequencies  while  the  disturbance 
bounds  dominate  in  the  high  frequencies.  However,  this  thesis  performs  a  digital 
QFT  design  and  the  sample  rate  places  certain  limits  on  the  ability  of  the  con- 
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troller  to  meet  specific  tracking  and  disturbance  bounds.  The  design  objectives 
are  therefore  redefined.  The  loop  transmission  function  is  designed  to  minimize 
the  phase  differential  between  the  loop  transmission  and  the  stability  bounds  while 
maximizing  the  gain  magnitude  in  the  low  frequencies.  Professor  Horowitz  refers 
to  this  approach  as  maximizing  the  amount  of  loop  transmission. 

While  performance  specifications  are  not  specifically  incorporated  into  the 
loop  shaping  process,  the  tracking  specifications  are  still  used  to  design  the  prefilter. 
The  tracking  specifications  for  a  fighter  type  aircraft  are  taken  from  previous  theses 
which  dealt  with  the  AFTI  F-16  and  are  in  appendix  D  [2,25  . 

Conclusion 

The  linear  equations  of  motion  describing  the  CRCA  across  a  range  of  sim¬ 
ulated  surface  failures  are  developed.  A  weighting  matrix  is  developed  which  con¬ 
verts  the  plant  transfer  function  matrix  (with  a  three-output  C  matrix)  into  a 
square  3x3  transfer  function  matrix. 


III.  The  Discrete  Plant  Model 


Introduction 

This  chapter  describes  the  transformations  of  the  plant  equations  required 
to  design  a  digital  controller  for  a  continuous  plant.  Certain  features  are  intro¬ 
duced  when  designing  a  control  system  for  a  sampled-data  system  which  are  not 
found  in  a  continuous  design.  These  unique  aspects  of  the  digital  approach  can 
be  traced  to  implementing  the  control  laws  as  difference  equations  in  a  computer 
and  to  introducing  a  sampler  and  zero  order  hold  (ZOH)  between  the  controller 
and  the  aircraft.  The  purpose  of  this  chapter  is  to  provide  a  brief  overview  of  the 
fundamentals  of  a  discrete  control  system,  examine  the  necessary  modifications  to 
applying  QFT  in  the  digital  domain,  and  the  development  of  the  discrete  CRCA 
model. 

Digital  System  Fundamentals 

The  z  Plane  A  digital  control  system  for  a  sampled-data  system  replaces  the 
analog  implementation  of  the  controllers  with  a  computer.  Fundamental  to  the 
use  of  a  computer  as  the  controller  for  a  continuous  dynamic  system  is  sampling. 
The  plant  outputs  must  be  sampled  to  convert  the  plant  response  to  a  digital 
representation.  At  the  same  time,  the  digital  output  of  the  computer  must  be 
converted  to  a  continuous  signal  which  is  then  applied  to  the  plant.  Both  processes 
affect  how  the  design  for  a  digital  controller  proceeds. 

The  sampling  process  converts  a  continuous  signal  into  an  impulse  train.  The 
value  of  the  pulse  at  discrete  intervals  is  given  by 

mn  =  -  kT"i  (3.1) 

where  Jd{kT)  is  the  discrete  value  of  the  continuous  function  e(f)  and  -  kT)  is 
the  unit  impulse  train.  The  Laplace  transform  of  the  expression  in  equation  (3.1) 
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is  given  by 

M*)  =  5£.«(HVr  (3-2) 

which  is  the  basis  for  the  transformation  which  maps  the  complex  variable  s  into 
the  complex  variable  z  by  the  relationship 

2  =  e‘T  (3.3) 

Such  a  transformation  is  called  the  Z  transform  and  is  a  mapping  of  a  discrete 
time  domain  function  into  the  complex  variable  z.  The  complex  variable  z  is  used 
in  a  representation  of  an  infinite  pulse  train  and  can  be  formulated  sometimes  in 
a  closed  form  or  always  in  an  open  form.  The  closed  form  can  be  manipulated 
as  a  transfer  function  which  lends  itself  to  design  work.  The  open  form  is  used 
to  implement  the  z  transfer  function  in  software  as  a  series  of  recursive  difference 
equations  [20,23], 

While  the  Z  transform  allows  the  tractable  manipulation  of  a  discrete  pulse 
train,  there  are  certain  aspects  of  the  z  plane  which  are  disadvantageous  to  the 
flight  control  system  designer  (24 j.  For  instance,  the  imaginary  axis  of  the  s  plane 
is  mapped  into  the  unit  circle  of  the  z  plane.  The  poles  of  the  plant  tend  to 
congregate  near  the  unit  circle  boundary  as  the  sampling  time  is  decreased.  The 
result  is  visibility  into  the  degree  of  stability  is  lost.  Also,  numerical  accuracy 
becomes  important  for  any  realistic  sampling  times;  the  difference  between  pole 
locations  becomes  extremely  small  as  sample  rates  decrease.  Therefore,  numerical 
accuracy  to  the  sixth  or  seventh  place  becomes  important.  In  addition,  the  f  and 
wn  contours  in  the  s  plane  become  quite  distorted  in  the  z  plane.  Again,  visibility 
into  system  performance  is  lost  when  the  interpretation  of  such  characteristics 
as  damping  or  undamped  natural  frequency  becomes  difficult.  Finally,  there  is 
no  clear  parallel  in  the  z  plane  for  such  classical  analysis  tools  as  root  locus  or 
frequency  response.  Fortunately  an  alternative  method  exists  which  allows  the 
analysis  of  a  sampled  data  system  which  does  not  have  these  limitations. 
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The  u'l  Plane  The  wt  plane  is  the  result  of  a  bilinear  transformation  which 
maps  the  complex  z  variable  into  the  complex  wt  variable.  Analogous  to  the 
complex  s  variable,  the  wt  variable  has  a  real  and  imaginary  part  given  by 

wt  =  u  -r  jv  (3-4) 

The  interior  of  the  unit  circle  is  mapped  into  the  left  side  of  the  familiar  cartesian 
coordinate  system.  The  exterior  of  the  unit  circle  is  mapped  into  the  right  side  of 
the  cartesian  coordinate  system.  The  wt  plane  and  the  z  plane  are  related  by 

T{z-\) 

(«>»-*)  '  ’ 

The  wt  plane  is  the  same  as  the  continuous  s  plane  in  form.  The  f  and  w„  contours 
once  again  have  the  familiar  shapes  and  techniques  such  as  Bode  analysis  are 
applicable  again.  In  other  words,  the  wt  plane  allows  the  control  designer  to 
analyze  a  sampled  data  system  using  the  continuous  domain  techniques.  However, 
there  is  no  free  lunch;  certain  features  appear  in  the  wt  plane  which  are  not  found 
in  the  s  plane. 

One  feature  of  the  wl  plane  not  found  in  the  -s  plane  is  the  dependance  of 
root  location  on  the  sampling  rate  of  the  digital  control  system.  As  described 
by  Professors  Houpis  and  Lamont  in  Digital  Control  Systems  [20],  the  correlation 
between  the  root  location  in  the  two  planes  is  dependent  on  sample  rates.  The 
imaginary  part  of  the  root  transforms  accurately  to  the  wt  plane  if  the  product 
of  the  s  plane  imaginary  part  and  half  the  sampling  rate  is  equal  to  or  less  than 
0.297.  Similarly,  if  the  square  of  the  product  of  the  s  plane  real  part  and  half  the 
sampling  rate  is  much  less  than  2,  the  real  part  of  the  root  shows  good  correlation 
in  both  planes.  If  either  of  the  above  conditions  is  not  true,  the  root  location  is 
warped  in  the  transformation  process  from  the  s  to  the  z  plane  to  the  wt  plane. 
The  degree  of  warping  can  be  further  increased  by  some  of  the  approximations  to 
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the  5  plane  to  z  plane  mapping.  If  the  warping  is  too  severe,  the  usefulness  of  the 
u't  plane  approach  to  digital  design  is  greatly  reduced  20  . 

A  second  feature  of  the  wi  plane  is  a  result  of  the  use  of  a  zero  order  hold 
(ZOH).  Continuous  dynamic  systems  such  as  aircraft  do  not  respond  well  when 
subjected  to  pulse  inputs  from  the  discrete  output  of  a  computer  or  sampler. 
Therefore,  a  ZOH  is  employed  to  convert  the  discrete  signals  into  a  continuous 
signal  20  .  The  largest  frequencies  present  in  the  reconstructed  signal  must  be 
roughly  eight  times  less  than  the  rate  at  which  discrete  signals  are  provided  to  the 
input  of  the  ZOH  in  order  to  minimize  jump  discontinuities  on  the  output  of  the 
ZOH  at  the  sample  instants.  The  ZOH  must  be  considered  part  of  the  basic  plant 
when  the  various  transformations  are  performed.  The  result  of  a  transformation 
from  the  sorz  plane  to  the  wi  plane  is  a  right-half-plane  zero  whose  value  is  given 
by 


A  minimum  phase  plant  in  the  sor  J  plane  is  therefore  converted  to  a  nonmini¬ 
mum  phase  plant  in  wl  plane.  While  a  nonminimum  phase  plant  is  not  inherently 
unstable  the  phase  response  of  the  plant  will  be  radically  different  from  a  similar 
minimum  phase  plant.  Consequently,  as  Professor  Horowitz  points  out,  the  ben¬ 
eficial  effects  of  feedback  are  reduced  and  the  designer  must  therefore  modify  the 
expected  performance  from  a  nonminimum  phase  plant. 

A  final  feature  of  the  w>  plane  is  due  the  characteristics  of  the  bilinear  trans¬ 
formation.  A  unequal  order  numerator  and  denominator  in  the  s  plane  is  trans¬ 
formed  through  the  z  plane  to  an  equal  order  numerator  and  denominator  in  the 
wt  plane.  For  instance,  if  a  s  plane  numerator  has  only  a  gain  term  and  the  de¬ 
nominator  has  four  poles  in  the  s  plane,  the  wi  plane  numerator  will  gain  four 
zeros:  at  least  one  right  half  plane  zero,  the  sampling  zero,  and  three  other  zeros, 
some  of  which  may  also  be  in  the  right  half  plane.  The  nonminimum  phase  zeros 
in  addition  to  the  sampling  zero  further  complicate  the  design  process  by  increas- 
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ing  the  severity  of  nonminimum  phase  characteristics.  An  equal  order  transfer 
function  also  has  frequency  characteristics  which  are  different  from  the  typical  s 
plane  transfer  function.  The  magnitude  of  the  wi  transfer  function  will  approach  a 
constant  value  rather  than  zero  as  the  frequency  increases.  At  the  same  time,  the 
phase  will  approach  zero  degrees  w’ith  increasing  frequency  (assuming  positive  gain 
and  minimum  phase).  The  frequencies  at  which  the  magnitude  becomes  constant 
and  the  phase  goes  to  zero  is  a  function  of  the  sampling  rate.  These  characteristics 
place  limits  on  the  performance  improvements  which  can  be  achieved  by  feedback. 
Therefore,  the  performance  limits  are  clearly  related  to  the  sampling  rate. 

Digital  QFT  The  introduction  of  the  nonminimum  phase  zeros  by  the  trans¬ 
formation  process  does  not  prevent  the  application  of  QFT  but  it  does  require  a 
certain  modification  of  the  loop  shaping  process.  Because  the  nominal  plant  is  non¬ 
minimum  phase,  the  nominal  loop  must  also  be  nonminimum  phase.  The  right-half 
plane  zeros  can  complicate  the  design  of  the  nominal  loop  transfer  function  because 
the  phase  response  characteristics  are  altered.  An  alternative  to  nonminimum 
phase  loopshaping  is  made  possible  by  the  introduction  of  an  all-pass-  filter  |2lj. 
A  function  is  initially  defined  as 

f  («■«)  =  (3.7) 

1  2 

which  is  just  another  expression  for  1.  The  loop  transfer  function  is  rearranged  so 
that  the  minimum  phase  part  of  the  loop  transfer  function  is  cascaded  with  the 
nonminimum  phase  sampling  zero. 

i™,  =  (£»,)  (l  -  ~)  (3  8) 

where  1  -  ^  is  the  sampling  zero. 

The  function  in  equation  (3.7)  is  cascaded  with  the  expression  in  equation  (3.8)  so 
that 

[Lnmrmw 0)  =  (L„„)  (l  -  ^)  (3.9) 
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Frequency 

Phase  Shift 

1  rad  sec 

0.72  degrees 

5  rad /sec 

3.58  degrees 

10  rad /sec 

7.21  degrees  ; 

25  rad  'sec 

17.8  degrees 

50  rad /sec 

34.7  degrees 

100  rad /sec 

64.0  degrees 

160  rad/sec 

90.0  degrees 

320  rad/sec 

127  degrees 

Table  3.1.  Phase  Shift  due  to  A(u>t) 


The  expression  is  rearranged  to 


L„mp  =  (L„r)  (l^)(L-g) 


(3.10) 


The  QFT  design  is  then  developed  based  on  the  gain  and  phase  characteristics 
of  the  minimum  phase  transfer  function,  Lmp ,  whose  relationship  to  the  orignal 
nonminimum  phase  transfer  function  is  given  by 


r  it  w,  Twi 

Lmp  =  [Lmp){  1  +  — -) 


(3.11) 


The  all-pass-filter,  A(u't),  is  given  by 


(3-12) 

I  2 

The  magnitude  of  A(ivl)  is  unity  over  the  frequency  range  but  the  phase  response 
is  given  by 


$  =  2  arc  tan 


(3.13) 


The  phase  shift  given  in  equation  (3.13)  is  no  longer  a  part  of  the  minimum  phase 
nominal  loop  transfer  function  but  mathematical  validity  requires  that  the  phase 
shift  be  accounted  for.  Instead,  the  phase  shift  shows  up  in  the  frequency  dependent 
stability  bounds.  Each  boundary  is  shifted  toward  the  zero  degree  phase  point  in 
accordance  with  the  expression  in  equation  (3.13).  For  the  sample  rate  used  in  this 
thesis,  the  phase  shift  due  to  the  all-pass-filter  is  given  in  Table  3.1.  Examination 
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of  Table  3.13  indicates  very  little  phase  shift  occurs  in  the  lower  frequencies.  Once 
the  phase  shift  becomes  significant,  the  boundaries  have  already  folded  into  the 
UHFB.  The  straight  boundary  of  the  UHFB  is  shifted  until  it  reaches  0  degrees. 

The  Hoffman  Algorithm  The  conversion  from  the  s  plane  to  the  wt  plane 
is  performed  using  the  Hoffman  algorithm  [9  .  The  algorithm  does  away  with  the 
intermediate  results  in  the  z  plane  and  therefore  avoids  the  numerical  problems 
involved  in  z  plane  analysis.  The  approach  solves  the  hyperbolic  tangent  function 
'  tanh(iy )  with  an  infinite  series  solution.  For  this  thesis,  the  series  was  trun¬ 
cated  at  the  fourth  term  jl  .  A  simulation  of  the  fourth  order  canard/flaperon 
actuators  is  performed  in  the  s  plane  and  the  wt  plane.  A  simulation  of  the  actua¬ 
tor  model  using  MatrixX  is  performed  first  with  a  continuous  step  input.  A  second 
simulation  with  a  continuous  step  input  is  performed  after  the  actuator  model  is 
transformed  to  the  wt  domain.  The  results  are  shown  in  Figure  3.1.  Examination 
of  the  results  shows  identical  responses  are  obtained  in  both  planes  indicating  the 
dominant  roots  are  not  warped  to  any  noticeable  extent. 

The  Transformed  Effective  Plant 

The  nominal  effective  plant  is  obtained  in  the  wt  plane  by  first  converting  the 
effective  MIMO  state  space  plant  while  still  in  the  s  plane  into  nine  SISO  transfer 
functions  using  the  expression  in  equation  (2.8).  The  actuator  transfer  function 
is  cascaded  with  each  effective  plant  transfer  function.  The  cascaded  transfer 
functions  are  then  converted  individually  back  into  the  state  space  format.  The 
Hoffman  algorithm  is  used  to  convert  the  nine  s  plane  state  space  models  into  nine 
wt  state  space  models.  Finally,  the  state  space  models  are  converted  back  into  nine 
transfer  functions  which  represent  the  nominal  effective  plant  in  the  wt  plane. 

Initially  the  fourth  order  actuators  are  cascaded  with  each  transfer  function 
before  transformation  to  the  wt  plane.  However,  nonminimum  phase  zeros  in  ad- 


Figure  3.1.  Step  Response  of  Fourth  Order  Actuator  in  s  and  wt  Plane 


dition  to  the  sampling  zero  are  introduced  by  the  conversion  process.  While  the 
nonminimum  phase  zeros  in  the  individual  elements  of  the  transfer  function  matrix 
are  not  intrinsically  detrimental,  the  right  half  plane  roots  of  the  determinant  of 
the  transfer  function  matrix  should  be  limited  to  the  sampling  zeros  in  order  to 
be  considered  minimum  phase.  In  this  case,  the  determinant  of  the  effective  plant 
had  nine  numerator  roots  consistently  appearing  in  the  right  half  plane:  three  at 
wt  -  160,  three  at  wt  =  289,  and  three  which  varied  around  wt  —  1100.  The  roots 
at  wt  =  289  and  wt  =  1100  are  attributed  to  the  fourth  order  actuators.  The 
all-pass-filter  could  be  modified  for  multiple  zeros  but  the  resultant  phase  shift  on 
the  stability  boundaries  would  become  too  severe.  Also,  the  roots  at  wt  =  1100 
varied  by  as  much  as  100  and  would  not  be  effectively  cancelled  by  the  all-pass-filter 
for  all  the  members  within  the  design  set.  Therefore,  the  fourth  order  actuators 
are  replaced  by  a  first  order  actuator  which  is  given  by 

(57.3)  •  (£-.)  (3.14) 

The  gain  term  premultiplying  the  transfer  function  is  necessary  to  convert  the 


output  of  the  model  from  radians  to  degrees.  The  control  input  matrix  is  in  terms 
of  degrees.  The  first  order  and  fourth  order  actuators  are  simulated  in  the  s  plane. 
The  results  are  shown  in  Figure  3.2.  Examination  of  the  results  shows  similar  rise 
times  but  the  first  order  response  lacks  the  time  lag  of  the  fourth  order  actuators. 


Figure  3.2.  Step  Response  of  Fourth  vs  First  Order  Actuator  in  s  Plane 

The  first  order  actuators  are  cascaded  with  each  transfer  function  in  the 
effective  plant  matrix  and  individually  transformed  into  the  wt  plane.  The  deter¬ 
minant  for  all  the  effective  plants  in  the  design  set  is  shown  to  be  minimum  phase 
excluding  the  nonminimum  phase  zeros  due  to  sampling  at  w  =  160  except  cases 
20,  22,  and  24.  These  three  cases  have  additional  zeros  and  are  therefore  excluded 
from  the  design  templates. 

The  Effective  Plant  in  the  wi  Plane  The  effective  transfer  function  plant 
matrix  in  the  wi  plane  for  the  ACM  entry  nominal  case  is  given  in  appendix  A. 


The  qXj  Terms  An  important  feature  of  a  QFT  MIMO  design  is  the  indi¬ 
vidual  design  of  each  feedback  loop.  In  order  to  accomplish  single  loop  design, 
the  nxn  MIMO  plant  is  mapped  into  n2  MISO  plant  equivalents.  The  following 
derivation  is  based  on  the  desire  to  formulate  a  simple  expression  for  the  desired 
compensator  [21  .  The  system  control  ratio,  T(w),  which  relates  the  vector  of 
inputs  to  the  vector  of  outputs  is  given  by 

T(tt>)  =  7  +  PG(w)  (3.15) 


The  control  ratio  for  any  single  input-output  pair  is  quite  complex.  However, if  the 
expression  is  rearranged  so  that 


P1  +  G(u>)]  T(w)  =  G(w) 


(3.16) 


then  only  the  inversion  of  the  plant  matrix  is  required.  Since  P(w)  is  a  matrix  of 
transfer  functions,  the  inverse  P(wt)  is  given  by 


pm  =  pH 

1  1  det\P(w)} 


(3.17) 


The  effective  plant  transfer  functions  are  formed  from  the  inverse  of  the  individual 
terms  of  inverse  of  the  plant  and  are  referred  to  as  the  qXJ  terms.  The  q,:  terms  are 
derived  by 

det\P(w)' 

%  =  <  3.18 

adj(p,}) 

While  simple  to  show,  performing  the  necessary  algebra  to  derive  the  qXJ  can  be 
tedious.  However,  the  nominal  plant  for  the  ACM  entry  condition  has  several 
fortuitously  placed  zero  terms  which  allows  a  general  formulation  of  the  q,j  terms 
to  be  derived.  Once  the  general  expression  is  derived,  the  individual  qXj  terms  are 
computed  from  the  coefficients  of  the  ptJ  terms.  The  derivation  is  started  with  the 
effective  plant  matrix. 

r  I 

a  0  0 

P*f!=  |  0  6  c  (3.19) 

Ode 


| 
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Where  a.  b.  c.  d.  and  e  are  the  individual  transfer  functions  of  the  effective  plant 
matrix.  The  inverse  of  a  square  matrix  is  given  by 

_  adj'P'fj 


PtS1  det\P'!f\ 
The  adjoint  of  Pe/f  is  derived  in  general  terms  by 


(3.20) 


adj  P tff  = 


0 


be  -  dc  0 

0  ae  —ac 

0  -ad  ab 


(3.21) 


The  determinant  of  Pt}j  is 


det'P'ff  j  =  a(6e  -  dc) 


(3.22) 


Therefore,  using  the  expression  in  equation  (3.18),  the  nominal  q,j  terms  are 


Q'if  = 


a 

oc 

oo 


oc 

oc 

bt-dc 

de-be 

t 

c 

dc-bt 

be-dc 

d 

b 

(3.23) 


The  nominal  Qejj  is  computed  and  the  results  are  found  in  the  appendix  E. 
Conclusion 

The  3x3  effective  plant  in  the  s  plane  is  transformed  to  the  wt  plane.  The 
transformed  plant  matrix  is  shown  to  be  minimum  phase  for  25  of  the  28  possible 
cases  in  the  design  set  for  ACM  entry  and  17  of  the  18  cases  for  TFTA.  The 
templates  which  are  used  to  determine  the  stability  bounds  are  based  on  a  reduced 
design  set  containing  the  25  minimum  phase  plant  transfer  function  matrices. 
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IV.  Compensator  and  Prefilter  Design 

Introduction 

The  three  compensators  and  three  prefilters  are  designed  using  QFT.  The 
first  loop  is  picked  on  the  basis  of  least  uncertainty  and  the  original  method  of 
QFT  is  used.  The  next  two  loops  are  designed  using  the  improved  method  of  QFT 
in  order  to  reduce  overdesign  of  the  compensators.  The  compensators  and  prefilters 
designed  in  the  wi  plane  are  transformed  to  the  z  plane  for  implementation  in  the 
aircraft  simulation. 

Design  of  Pitch  Rate  Loop 

Compensator  One  The  first  loop  to  be  shaped  is  the  pitch  rate  feedback  loop. 
The  decision  which  loop  to  design  first  is  important  because  that  loop  compensator 
will  exhibit  the  most  overdesign.  After  the  first  loop  compensator  is  designed,  that 
compensator  is  incorporated  into  the  nominal  loop  transmission  function  of  the 
next  two  loops  when  the  improved  method  of  QFT  is  applied.  Therefore,  the 
decision  which  loop  to  shape  first  colors  the  remainder  of  the  design.  However,  no 
consensus  appears  to  exist  on  how  to  pick  the  first  loop.  Previous  thesis  efforts 
have  started  with  both  the  most  uncertain  and  the  least  uncertain  transfer  function 
set.  Neither  approach  has  been  shown  to  be  the  optimal  approach.  For  this  thesis, 
based  on  advice  from  Professor  Horowitz,  the  least  uncertain  transfer  function  set 
is  picked  as  the  first  loop.  The  level  of  uncertainty  is  indicated  by  the  size  of  the 
templates  as  determined  by  the  phase  angle  and  magnitude  data  for  q,t.  In  this 
case,  the  smallest  templates  and  therefore  least  uncertain  are  obtained  for  qu- 

The  compensator  is  developed  by  adjusting  the  frequency  response  of  the 
loop  transmission  function  (loop  shaping)  as  depicted  on  the  Nichols  chart.  The 
loop  frequency  response  is  adjusted  by  locating  additional  poles  and  zeros  which 
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are  cascaded  with  the  unimpaired  case.  gi}(ir/).  so  that  the  frequency  response  of 
the  loop  at  a  particular  frequency  has  the  desired  phase  and  gain  characteristics. 
The  process  of  loop  shaping  facilitates  the  point-wise  tailoring  of  the  frequency 
response  of  the  compensator  and  plant  as  a  whole  unit.  The  transfer  function  used 
to  start  the  loop  shaping  is  called  the  nominal  loop  transmission  nonminimum 
phase  transfer  function  and  is  given  by  L\{wt)nmp  =  gn(u>f)nmp.  For  the  first  loop, 
this  function  is 

_  1.4981  •  10~3wf(wt  +  G.00969)(tf/  -  1.07)(u>/  -  3723 )(wi  -  160) 

iltu'Jnmp  -  —  —  —  ^  3.370)(U>/  -  1.649)(iet  ~  0.00402  =  j0.03561) 

(4.1) 

which  is  in  the  form  of  equation  3.8.  The  nominal  loop  transmission  function  is 
both  unstable  and  nonminimum  phase.  Using  an  all-pass-filter,  the  nonminimum 
phase  nominal  transfer  function  is  transformed  to  the  minimum  phase  nominal 
case  which  is  in  the  form  of  equation  3.10.  The  resulting  transfer  function  used  to 
begin  the  loop  shaping  process,  L\MP  =  qn ,  is  given  by 

1. 4981-  10“ 3wi(wi-  0.00969) (art  +  1.07) (u;t  +  3723) (u>»  +  160) 
lMP  ~  (wt  -  19.9) (wl  t  3.370) (u»f  -  1.649) [wi  -  0.00402  ±  jO.03561)  1  ' 2 ' 

Once  a  loop  is  shaped  which  meets  the  requirements,  the  compensator,Gi(u^t), 
is  obtained  by  multiplying  the  loop  transmission  function  by  the  inverse  of  the 
nominal  transfer  function. 

The  loop  shaping  is  performed  based  on  the  stability  bounds  as  determined 
by  the  templates  for  gjj(u>/).  The  constant  magnitude  circle  of  3  dB  on  the  Nichols 
chart  is  picked  as  the  absolute  stability  region  which  the  loop  transmission  function 
can  never  penetrate.  The  3  dB  limit  corresponds  to  a  closed  loop  peak  magnitude 
response  for  a  step  input  of  1.41.  The  templates  are  translated  but  never  rotated 
about  the  absolute  stability  region.  The  nominal  plant  is  marked  on  the  chart 
and  the  marks  are  connected  to  form  the  boundary  for  a  particular  frequency.  At 
the  low  frequencies,  the  templates  are  sufficient  to  trace  the  boundaries  because 
negligible  phase  shift  is  introduced  by  the  all-pass-filter.  At  frequencies  above  6 
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rad/sec,  the  phase  shift  becomes  significant  but  at  the  same  time,  the  templates 
have  shrunk  to  a  constant  phase  with  a  length  equal  to  12  dB.  The  stability  bound¬ 
aries  therefore  become  overlapping  vertical  lines  tangent  to  the  constant  magnitude 
circle  in  the  minimum  phase  case.  The  resulting  single  boundary  is  called  the  uni¬ 
versal  high  frequency  boundary  (UHFB).  The  UHFB  is  shifted  toward  the  zero 
degree  line  on  the  Nichols  chart  with  increasing  frequency  due  to  the  effect  of  the 
all-pass-filter. 

The  loop  transmission  function  must  continually  be  to  the  right  of  the  stabil¬ 
ity  bound  at  each  frequency  and  must  be  below  the  UHFB  before  the  zero  degree 
phase  point  is  reached.  The  requirement  to  be  below  the  UHFB  at  0  degrees  places 
a  limit  on  the  amount  of  loop  transmission  which  can  be  achieved  in  the  design. 
For  the  first  loop  of  this  design,  the  task  is  eased  somewhat  because  the  nomi¬ 
nal  point  on  the  templates  migrated  to  the  top  of  the  templates  with  increasing 
frequency.  Therefore,  the  UHFB  is  above  the  constant  magnitude  circle  and  thus 
greatly  eased  the  task  of  avoiding  the  UHFB  in  the  higher  frequencies. 

The  objective  of  loop  shaping  in  a  digital  design  is  to  maximize  loop  trans¬ 
mission  magnitude  while  satisfying  the  boundary  constraints  at  each  frequency. 
Boundary  constraints  are  satisfied  by  minimizing  the  phase  difference  between  the 
loop  transmission  function  and  the  stability  boundaries.  This  approach  maximizes 
system  performance  within  the  constraints  due  to  sampling.  The  loop  shaping  is 
performed  with  the  minimal  number  of  poles  and  zeros  possible  in  order  to  simplify 
the  resulting  compensator.  Since  this  design  is  performed  in  the  w>  plane,  an  equal 
number  of  poles  and  zeros  must  be  added  to  L\MP{w).  The  loop  transmission  gain 
is  increased  until  the  stability  boundaries  are  just  met  in  the  higher  frequencies. 
The  result  of  maximizing  the  loop  transmission  is  a  high  crossover  point  in  the  wi 
plane.  The  loop  transmission  function  is  plotted  in  figure  4.1  and  the  data  used  to 
plot  the  loop  transmission  is  found  in  appendix  F. 


Since  the  loop  transmission  function  is  equal  to  the  minimum  phase  nominal 
plant  multiplied  by  the  compensator,  division  is  performed  to  obtain  the  compen¬ 
sator.  As  an  aside,  the  resulting  compensator  is  only  applicable  to  the  minimum 
phase  plant;  however,  the  actual  plant  to  be  controlled  is  nonminimum  phase.  A 
polar  plot  analysis  reveals  a  negative  gain  is  required  to  maintain  the  necessary 
loop  transmission  characteristics  for  the  nonminimum  phase  plant.  Thus,  for 

L\ MP[wi)  =  9ii(u»0  •  Gi(wi)  (4.3) 


After  multiplication  by  a  negative  one  to  account  for  the  effects  of  the  all-pass-filter, 
the  compensator  for  the  nonminimum  phase  pitch  rate  loop  is 


e  ,(«>*  + 0.005)(u>f  +  4.6) {wl  +  44.5) (wi  +  12000)  „  c, 

Gl{w,)  ~  “215-4 - («7)(^f-  -TT.7)(^f  “SQOOK^r  ^5047) -  (4’5) 


Pre filter  One  The  purpose  of  the  prefilter  is  to  position  the  frequency  re¬ 
sponse  of  the  loop  transmission  function  within  the  frequency  bounds  as  described 
by  the  upper  and  lower  tracking  bounds.  Since  the  compensator  is  designed  to 
maximize  the  loop  transmission,  the  closed  loop  time  response  to  a  step  input 
without  a  prefilter  is  unrealistically  fast.  The  prefilter  is  designed  as  a  low  pass  fil¬ 
ter  which  removes  the  higher  order  frequencies  of  the  pilot  command  and  therefore 
slows  down  the  resulting  closed  loop  system  response.  At  the  same  time,  the  closed 
loop  frequency  response  must  remain  within  the  bounds  for  the  entire  design  set. 
Therefore,  the  worst  case  frequency  responses  for  each  impairment  in  the  design 
set  is  used  to  form  an  upper  and  lower  response  boundary.  The  desired  upper  and 
lower  tracking  bounds  are  given  in  the  form  of  transfer  functions  and  are  found  in 
the  appendix  A.  The  prefilter  is  then  used  to  position  the  closed  loop  responses 
of  all  the  members  of  the  design  set  within  the  desired  response. 
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Deriving  the  frequency  response  of  the  closed  loop  transfer  functions  for  the 
25  transfer  function  of  the  design  set  is  somewhat  tedious  so  numerical  methods 
are  used  instead.  In  general,  the  closed  loop  transfer  function  for  a  system  with 
unity  feedback  without  a  prefilter  is  given  by 


T,(wi)  = 


G(u'/)q„(w>) 


1  —  G(wf)q,(wf) 

where  T,(wf)  and  q,(wf)  are  the  closed  loop  transfer  functions  and  corresponding 
open  loop  transfer  functions,  respectively,  for  the  design  set  and  G(wi)  is  the  com¬ 
pensator.  The  magnitude  of  the  frequency  response  for  each  closed  loop  transfer 
function  of  the  design  set  is  found  in  decibels  in  terms  of  the  open  loop  compensator 
and  each  plant  in  the  design  set  as  follows 

T.  -  20  log  [  » W  ,1  (4.7) 

[1  -  G(jv)q,{jv)\ 

T,  =  20 log  'G(jv)  ■+■  20 log  q,(jv)  -  20 log  "'1  -  G(jv)q,(jv)  '  (4.8) 


The  last  equation  is  solved  repeatedly  in  a  loop  for  each  plant  case  and  the  max¬ 
imum  and  minimum  bounds  for  a  discrete  set  of  frequencies  in  the  wt  plane  are 
determined.  The  maximum  upper  case  response  is  subtracted  from  the  desired 
upper  tracking  bounds  and  the  same  operation  is  performed  on  the  lower  set  of 
bounds.  The  result  is  a  new  set  of  bounds  on  F(jv)  which  are  used  to  develop 
the  prefilter.  The  prefilter  guarantees  that  the  response  of  the  entire  design  set  is 
constrained  within  the  desired  tracking  bounds.  Therefore,  the  magnitude  of  the 
closed  loop  transfer  function,  T(jv),  cascaded  with  the  is  the  prefilter,  F(jv),  is 
given  by 

\Ti{jv)  (=  |r(jt;)||F(yt,)'  (4.9) 


(4.10) 


Using  the  above  approach,  the  prefilter  for  the  first  loop  is  given  by 

_  0.0005156(^400)^-160) 
v  ’  Wi  +  3)(urt  +  11) 

Figure  4.2  displays  the  tracking  response  of  the  nominal  closed  loop  response  to 
a  step  input  in  the  wt  plane.  The  nominal  closed  loop  system  consists  of  the 
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compensator  and  the  9n„mr(u'f)  with  unity  feedback  and  the  prefilter  on  the  input. 
This  single  input,  single  output  simulation  is  performed  as  a  quick  check  on  both 
stability  and  tracking  within  the  desired  tracking  boundaries  and  is  not  intended 
to  replace  a  full-up  simulation.  As  can  be  seen,  the  step  response  is  within  the 
desired  upper  and  lower  boundaries. 


Figure  4.2.  Closed  Loop  Step  Response  With  Prefilter  for  Pitch  Rate  Loop 


Design  of  Roll  Rate  Loop 

Compensator  Two  The  second  loop  designed  is  the  roll  rate  loop.  The  tem¬ 
plates  are  generated  using  both  the  original  method  and  the  improved  method. 
The  original  method  develops  the  templates  based  on  the  frequency  response  of 
the  922  transfer  function.  The  improved  method  develops  the  templates  by  in¬ 
corporating  information  from  the  first  loop  design  into  the  second  loop  using  the 
following  equations 

L2tn  (u>t)nmp  =  92922,  [wl)nnp  (4.1l) 

L2.,{wt)nmp  =  g2?22  (wOnmp-.-^t  ~  (4.12) 

1  +  l>i  ~  7 


4-7 


where 


and 


q22.{W,)nmp  =  922  («’') 


nmp  . 


i-Ii 


1  ~ 


<hlQ22 

912921 


(4.13) 


(4.14) 


Comparison  of  the  template  data  generated  by  q22  and  922,  reveals  some  changes  in 
template  sizes  in  the  lower  frequencies  but  identical  results  in  the  higher  frequen¬ 
cies.  The  conclusion  drawn  is  since  the  linear  aircraft  model  is  inertially  decoupled, 
there  is  little  to  be  gained  by  using  the  improved  method  between  the  longitudinal 
pitch  mode  and  lateral  roll  mode.  However,  the  improved  method  templates  are 
still  used  in  this  thesis  for  two  reasons.  First,  the  control  surfaces  are  not  decou¬ 
pled;  eight  of  the  nine  surfaces  provide  inputs  to  both  the  longitudinal  and  lateral 
modes.  Second,  inertial  coupling  is  observed  to  occur  for  certain  control  surface 
failures. 


The  original  922(w/),,mp  transfer  function  is  used  as  the  nominal  loop  nonmin¬ 
imum  phase  transmission  loop  because  'y  is  negligible  throughout  the  frequency 
range  for  the  nominal  case.  The  g22nmF  transfer  function  is  given  by 

1.0603  •  10~3wi(wt  r  3489)  (wt  -  160)  -  0.0628  ±  j0.007) 

2(wi)nmp  -  (u„  _  I9.9)(u;/  j  1.7l3)(u>/  -  0.0562)(wr  +  0.1785  ±  j0.8124)  ^  '  ' 

The  nominal  transfer  function  is  unstable  and  nonminimum  phase.  The  all-pass-filter 
method  is  applied  to  nonminimum  phase  loop  transmission.  The  resulting  mini¬ 
mum  phase  transfer  funciion,  q22  =  L2mp  used  in  the  loop  shaping  process  is  given 
by 

1.0603  •  10~3wt(wf  +  3489)(wt  +  160)(wf  +  0.0628  ±  ,70.007) 

2mp  ~  (wt^  19. 9) (wt  +  1.7 13) (wt  -  0.0562) (wt  +  0.1785  ±  j0.8124)  '  ' 

The  3  dB  constant  magnitude  circle  is  again  used  to  define  the  absolute  stability 
bounds.  As  in  the  first  loop  design,  the  roll  rate  loop  is  shaped  to  maximize 
low  frequency  gain  within  the  constraints  of  meeting  the  high  frequency  stability 
boundaries  in  the  wt  plane.  The  resultant  high  crossover  value  is  similar  to  the  first 
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loop  transmission  crossover.  Dividing  the  final  minimum  phase  loop  transmission 
by  the  nominal  loop,  L2MP(w),  and  changing  the  sign  on  the  gain  term  to  account 
for  the  all-pass-filter  yields  the  following  compensator  for  the  nonminimum  phase 
roll  rate  loop. 


G2(u-f)  =  -252.1 


(wl  +  0.8) (wl  -+-  7.09) (wt  45.9) (wl  -t-  749) (wl  14000) 
(wl)(wi  i-  0.305) (wl  -r  1440) (wl  +  2500) (wl  -r  8000) 


(4.17) 


Pre filter  Two  Applying  the  same  method  used  to  design  the  first  prefilter, 
the  second  set  of  prefilter  bounds  are  determined.  The  resulting  prefilter  is  given 
by 

0.0001  (wf  +  500 )(u»  -  160) 

(wi  —  2)(u>f  —  4) 

The  closed  loop  step  response  for  the  compensator  and  q22„nr,(u,,)  is  simulated  sim¬ 
ilarly  to  the  first  loop  design.  The  tracking  step  response  is  displayed  in  Figure  4.3 
and  is  seen  to  be  within  the  desired  tracking  upper  and  lower  bounds. 


F2(u’/) 


Figure  4.3.  Closed  Loop  Step  Response  With  Prefilter  for  Roll  Rate  Loop 


Design  of  Yaw  Rate  Loop 


Compensator  Three  The  improved  method  of  QFT  is  used  to  design  the 
compensator  for  the  last  channel  which  controls  the  yaw  rate  variable.  While 
there  is  no  inertial  coupling  between  the  longitudinal  and  lateral  modes  for  the 
unimpaired  case,  significant  coupling  exists  between  the  two  lateral  modes.  The 
inertial  coupling  is  due  to  the  off  diagonal  terms.  923(H  and  931  (H-  of  the  transfer 
function  matrix.  The  templates  are  developed  by  modifying  the  933(u>f)nmp  term 
to  incorporate  the  all-pass-filter  and  previous  loop  compensators  into  the  transfer 
function.  The  expression  for  the  effective  q3z,(wi)nmp  is  given  by 


933,  (Hr 


933  ( 

f  -  A 


(4.19) 


where 


9.165-  10-4w/{wi  +  3757) (wf  -  160) (wf  +  0.0628  ±>0.007) 
3(H„mp(H  -  (u.,  19.9) (wf  -  0.00125)  (wf  -  0.0562)  (wf  -r  0.2763  ±  1.89) 

(4.20) 


?  —  (1  +  -^i)(l  +  L2)  -  I12 
A  =  ">23(1  —  Li)  —  113(1  ^  L2 )  —  (112U2  +  I13U3) 


(4.21) 

(4.22) 


However,  due  to  the  decoupling  of  the  longitudinal  and  lateral  modes,  i12  and  i13 
are  zero.  After  some  algebraic  manipulation,  the  expression  for  933,  is  reduced  to 


933,  (*^)r»mp  — 


933(Hnmp(l  ^  H 
1  -t-  L2  -  123 


where 


I23  = 


(922K933) 


(4.23) 


(4.24) 


(923)  (932 ) 

After  the  algebraic  manipulations  are  complete,  the  effective  nonminimum  phase 
nominal  plant  transfer  function  for  933  is  given  by 


933,(H"mp  =  933  (Hr 


(wf  +  0.06268  ±  >0.00992) 
(wf  +  0.0628  ±  >0.007) 


(4.25) 


1 


,  A  9.165  •  10_4tr/(u-»  -  3757) (wt  -  160) (it-/  -  0.06268  ±  j0.0992) 

933 .K’Onmp  -  ^  —  19  _  o.00125)(wf  -  0.0562)H  Z 0.2763  ±  1.89) 

(4.26) 


After  modification  by  the  all-pass-filter,  the  resulting  nominal  minimum 
phase  loop  transmission,  LzMP,  is  used  to  develop  the  templates  and  start  the 
final  loop  shaping.  The  loop  is  shaped  as  in  the  first  two  loops  and  the  resulting 
compensator  for  the  yaw  rate  channel  is  given  by 

G  .  v  _26  >''-0.l)(u’^  2.45)(u’>  -i-  3.55)(u;t  30) 

3  (u>f )(«.’/  ■*  350) (wt  -  0.025  —  jO.0968) 


(4.27) 


Prefilter  Three  Using  the  same  upper  and  lower  tracking  response  bounds 
as  used  for  the  roll  rate  channel,  the  prefilter  for  the  third  loop  is  developed  using 
the  same  procedures  as  previously  employed.  The  third  prefilter  is  given  by 


■F3(teO  = 


0.001  (u»  -  500 )(wi  -  160) 
(wt  +  2)(tet  +  4) 


(4.28) 


A  single  input  single  output  simulation  of  the  closed  loop  system  is  performed  as 
in  the  previous  two  loops.  The  resulting  step  tracking  response  of  the  compensator 
with  the  nominial  q33nfnp(wt)  and  the  prefilter  is  shown  in  figure  4.4. 


Transformation  to  the  z  Plane 


While  the  design  is  performed  in  the  wt  plane,  the  compensators  and  pre- 
filters  are  implemented  in  the  aircraft  simulation  in  the  z  plane.  The  bilinear 
transformation  is  used  to  convert  the  transfer  functions  into  the  z  plane.  The 
relationship  between  the  z  and  the  wt  planes  is  given  by 


wt  = 


2  z-  1 

fz+1 


(4.29) 


where  T  is  the  sampling  time  of  0.0125  sec. 
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by 


Figure  4.4.  Closed  Loop  Step  Response  With  Prefilter  for  Yaw  Rate  Loop 
Using  the  bilinear  conversion,  the  compensator  for  the  pitch  rate  loop  is  given 


,  _  —20.2983(2  -  0.564792) (2  4-  0.973684)(2  -  0.944107)(2  -  .999938) 

“  (2  -  0.978973 )(2  -  0.937984) (2  +  0.938544)(2  -  1) 

(4.30) 


The  prefilter  for  the  pitch  rate  loop  is  given  by 

0.003315(2  -1-  0.428571) 


*i (*)  = 


(4.31) 


(2  -  0.871345)(2  -  0.963190) 

The  compensator  for  the  roll  rate  loop  is  given  by 

_  -20.1537(2  -  0.554153)(2  +  0.647963)(2  +  0.977401)(2  -  0.915136)(2  -  0.990050) 
2'Z’  _  (2  -  l) (2  4-  0.800000)(z  +  0.879699)(2  +  0.960784)(2  -  0.996195) 


and  the  prefilter  for  the  loop  is  given  by 

0.000794941(2  +  0.515152) 


F2(z)  = 


(2  -  0.951220)(z  -  0.975309) 

The  compensator  for  the  yaw  rate  loop  is  given  by 

-10.8538(2  -  0.684211)(z  -  0.969837)(2  -  0.956588)(2  -  0.998751) 


(4.32) 

(4.33) 


c»W  = 


(2  -  1)(2  +  0.372549)(2  -  0.999687  ±  j0. 00120962) 


(4.34) 
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and  the  prefilter  is 


(4.35) 


0.000794941(2  -  0.515152) 
~  (z  -  0.951220) (2  -  0.975309) 


Conclusion 


Three  rate  controllers  and  associated  prefilters  are  designed  in  the  wt  plane. 
Each  controller  and  prefilter  is  simulated  with  the  nominal  q„  and  is  shown  to 
have  no  overshoot  and  zero  steady  state  error.  The  controllers  and  prefilters  are 
transformed  from  the  u>t  plane  to  the  z  plane  for  implementation  in  the  simulation. 
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V.  Aircraft  Simulation  Setup  and  Results 


Introduction 

The  QFT-designed  digital  FCS  for  the  CRCA  is  simulated  with  the  three  pilot 
rate  commands.  Single  and  multiple  pilot  commands  are  simulated.  The  closed 
loop  system  is  shown  to  be  robust  for  the  complete  range  of  failures  modeled  by 
the  reduced  design  set  for  ACM  entry.  Linearity  of  response  is  shown  for  a  large 
range  in  magnitude  of  the  command  input.  The  TFTA  condition  is  simulated 
with  the  pitch  rate  command  but  is  found  to  be  unsatisfactory.  The  first  order 
actuator  model  is  replaced  with  the  fourth  order  model.  Identical  digital  filters 
are  cascaded  with  the  digital  controllers  to  account  for  the  change  in  phase  angle 
response  of  the  loop  transmission  due  to  the  additional  poles  of  the  fourth  order 
actuators.  The  closed  loop  response  is  seen  to  be  stable  and  almost  equals  the 
performance  obtained  with  first  order  actuators  in  the  longitudinal  mode.  The 
lateral  mode  exhibits  a  greater  oscillatory  behavior  than  previously  obtained  with 
the  first  order  actuators.  The  output  of  the  weighting  matrix  is  examined  for  all 
nine  control  surfaces.  The  effects  of  saturation  on  the  continued  control  surface 
cooperation  are  studied. 

Simulation  Setup 

The  CRCA  is  simulated  using  the  nonlinear  dynamic  simulation  capabilities 
of  MatrixX.  While  the  CRCA  is  simulated  using  the  linear  perturbation  equations 
employed  in  the  design  process,  the  actuators  are  modeled  to  include  deflection 
and  rate  limits  which  are  not  part  of  the  design  equations.  The  closed  loop  system 
is  simulated  as  a  sampled  data  system.  The  aircraft  and  actuators  are  modeled 
as  continuous  systems  while  the  prefilters  and  controllers  are  modeled  as  discrete 
systems  with  a  sample  rate  of  80  Hz.  The  sampler  and  ZOH  are  provided  as  part  of 
the  discrete  block  in  the  library  of  MatrixX  simulation  options.  The  block  diagram 
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Figure  5.1.  Simulation  Setup 

of  the  closed  loop  system  is  shown  in  Figure  5.1.  Commands  in  the  form  of  pulses 
are  provided  to  the  three  inputs.  The  weighting  matrix  distributes  the  commands 
inputs  to  the  nine  actuators.  All  eight  states,  control  deflections,  control  rates, 
and  the  generated  states  of  longitudinal  and  lateral  accelerations  are  provided  as 
output. 

The  simulations  for  each  command  input  are  performed  in  a  loop.  The 
nominal  aircraft  is  simulated  and  the  resulting  states,  actuator  deflections,  and 
actuator  rates  are  stored  in  a  matrix.  The  nominal  aircraft  is  replaced  by  an  aircraft 
with  simulated  control  surface  damage  (case  2  in  appendix  B)  and  the  simulation 
is  repeated.  Similarly,  the  rest  of  the  failure  cases  in  the  reduced  design  set  for 
the  ACM  entry  condition  are  sequentially  simulated.  The  nominal  CRCA  for  the 
TFTA  condition  is  simulated  but  is  found  to  be  unsatisfactory.  The  unsatisfactory 
attempt  to  extend  the  design  to  the  TFTA  condition  is  discussed  further  in  the 
next  chapter. 

Interpretation  of  the  simulation  results  must  be  done  while  remembering  that 
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the  linear  equations  of  motion  are  developed  using  perturbation  methods.  The 
linear  equations  ofl  motion  are  developed  with  the  assumption  that  the  approxi¬ 
mations  sinx  =  x  or  cosx  =  1  are  valid.  If  small  angle  assumptions  are  exceeded, 
the  simulation  results  are  not  valid.  Therefore,  the  initial  transient  state  responses 
can  be  expected  to  closely  reflect  reality  within  the  limits  of  the  equations  of  mo¬ 
tion.  On  the  other  hand,  if  a  new  trim  condition  is  established  at  the  completion 
of  a  command,  the  system  response  can  not  be  expected  to  closely  approximate 
reality.  For  instance,  a  simulation  involving  a  pulse  roll  command  leaves  the  air¬ 
craft  in  a  banked  orientation  since  the  aircraft  does  not  have  a  restoring  moment 
which  would  otherwise  reestablish  the  initial  trim  conditions.  On  the  other  hand, 
the  aircraft  tends  to  return  to  the  initial  trim  conditions  after  the  completion  of  a 
pulse  pitch  command  and  therefore  the  solutions  to  the  equations  of  motion  tend 
to  remaim  within  the  confines  of  the  initial  assumptions. 

Actuators  The  control  surface  actuators  are  modeled  using  a  deflection  lim¬ 
iter  and  a  limited  integrator.  In  addition  to  acting  as  a  rate  limiter,  the  limited  in¬ 
tegrator  also  discontinues  integrating  whenever  the  rate  limits  are  exceeded.  When 
the  actuator  rates  are  no  longer  exceeded,  the  model  immediately  starts  integrat¬ 
ing  again  and  therefore  prevents  a  lag  due  to  “integrator  windup”  in  the  actuator 
response.  The  values  for  the  deflection  and  rate  limits  are  given  in  Table  2.1. 

Commands  The  simulated  pilot  commands  are  pitch  rate,  roll  rate,  and  yaw 
rate.  The  commands  are  ramped  to  avoid  saturating  the  actuator  models.  The 
commands  reach  full  value  or  return  to  zero  in  0.2  seconds. 

Single  Command  Results 

Each  command  is  simulated  twice  with  the  reduced  design  set  for  the  ACM 
entry  condition.  The  simulation  is  performed  first  with  a  small  magnitude  com¬ 
mand.  The  objective  is  to  study  the  aircraft  performance  within  a  linear  region 
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of  response.  The  second  simulation  is  performed  using  a  large  magnitude  com¬ 
mand.  The  purpose  of  simulating  a  command  which  saturates  or  nearly  saturates 
the  actuator  models  is  to  study  the  nonlinear  response  of  the  linearly  designed 
control  system.  The  left  canard,  left  outer  flap,  left  inner  flap,  left  elevator,  and 
rudder  actuator  deflection  and  rate  responses  are  plotted  once  in  the  single  input 
simulations  to  demonstrate  typical  deflections  and  velocities.  The  states  which 
are  plotted  for  each  simulation  run  are  pitch  rate,  roll  rate,  yaw  rate,  pitch  angle, 
bank  angle,  sideslip,  normal  acceleration,  and  lateral  acceleration.  The  results  of 
simulating  all  the  cases  within  the  reduced  design  set  are  found  in  appendix  G. 

Pitch  Rate  Command  for  ACM  Entry  The  first  pitch  rate  command  sim¬ 
ulated  is  a  2  degree/second  pitch  command  for  three  seconds.  The  command 
generates  1.96  g  at  the  center  of  gravity  (eg)  of  the  aircraft.  The  response  plots  for 
all  25  cases  can  be  found  in  Figure  G.l  through  Figure  G.13  of  appendix  G.  The 
surface  deflections  and  rates  are  well  within  their  limits.  After  the  command,  the 
surfaces  are  seen  to  return  to  their  trimmed  positions.  The  pitch  rate  state  exhibits 
no  overshoot  and  is  seen  to  be  identical  for  all  aircraft  in  the  reduced  design  set. 
The  commanded  pitch  angle  also  exhibits  no  overshoot  and  is  identical  for  all  25 
cases  in  the  reduced  design  set.  The  uncommanded  rates  are  shown  to  be  nearly 
zero  and  the  associated  angles  are  also  virtually  zero.  The  normal  acceleration  is 
shown  to  increase  smoothly  with  no  overshoot  and  builds  to  a  full  1.9  g.  Lateral 
acceleration  is  quite  small.  The  results  for  simulating  the  unimpaired  aircraft  with 
the  pitch  rate  command  is  shown  in  figures  5.2  and  5.3. 

The  pitch  rate  command  is  simulated  for  the  second  time  with  a  magnitude  of 
12  degrees/second  for  three  seconds.  The  command  generates  6.8  g  at  the  eg  of  the 
aircraft.  The  response  plots  for  all  25  cases  can  be  found  in  Figure  G.14  through 
Figure  G.21  of  appendix  G.  No  saturation  is  seen  in  any  of  the  surfaces.  The 
peak  rate  required  to  perform  the  command  is  55  degrees/second  and  the  largest 
deflection  is  13  degrees.  The  commanded  pitch  rate  state  and  associated  pitch 


Figure  5.2.  Pitch,  Roll,  and  Yaw  Rate  Response  for  2  DPS  Pitch  Rate  Input 


angle  are  identical  for  all  25  aircraft  and  no  overshoot  occurs.  The  uncommanded 
rates  are  shown  to  be  virtually  decoupled.  Two  damage  cases  are  seen  to  maintain 
about  0.5  degree  of  sideslip  after  the  initial  transient  while  the  remainder  of  the 
damage  cases  show  virtually  no  sideslip.  No  bank  angle  is  maintained  in  any  of  the 
cases  after  the  transients  die  out.  Normal  acceleration  increases  smoothly  but  does 
not  reach  the  full  6.8  g  commanded.  Negligible  lateral  acceleration  is  observed. 

Roll  Rate.  Command  for  A  CM  Entry  The  roll  rate  command  is  initially  sim¬ 
ulated  with  a  35  degree/second  magnitude  for  0.5  seconds.  Surface  position  deflec¬ 
tions  and  rates  are  seen  to  be  well  within  limits  except  for  three  cases:  6,7  and  25. 
The  exception  cases  demonstrate  much  larger  surface  movements  and  saturation 
is  seen  in  the  rudder  for  cases  6  and  7.  The  surfaces  continue  to  deflect  after  the 
command  is  finished.  Examination  of  the  uncommanded  states  indicates  a  larger 
drift  from  the  equilibrium  position  for  cases  6  and  7  than  for  the  rest  of  the  cases. 
The  commanded  roll  rate  and  bank  angle  are  identical  for  all  the  cases  simulated. 
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Figure  5.3.  Pitch  Angle,  Bank  Angle,  and  Sideslip  for  2  DPS  Pitch  Rate  Input 

No  overshoot  is  observed  in  the  roll  rate  state  or  bank  angle.  Cross-coupling  to 
the  longitudinal  pitch  rate  state  is  limited  to  less  than  0.08  degrees /second  for  24 
of  25  cases.  Coupling  to  the  lateral  yaw  rate  is  limited  to  some  transient  behavior 
before  being  driven  to  zero.  The  response  plots  for  all  25  cases  can  be  found  in 
Figure  G.22  through  Figure  G.29  of  appendix  G. 

The  second  roll  rate  command  simulated  is  135  degree  'second  for  0.5  seconds. 
Extensive  saturation  is  seen  in  all  the  control  surfaces  after  the  pulse  command  is 
complete.  Cases  6  and  7  are  seen  to  diverge  completely  due  to  loss  of  control  in 
the  uncommanded  states.  The  commanded  roll  rate  and  bank  angle  are  identical 
for  all  cases  for  the  first  10  seconds  of  the  simulation.  No  overshoot  is  observed 
in  pitch  rate  or  bank  angle.  Cross-coupling  is  increased  due  to  the  loss  of  con¬ 
trol  of  uncommanded  states.  The  response  plots  for  all  25  cases  can  be  found  in 
Figure  G.30  through  Figure  G.37  of  appendix  G. 
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Yaw  Rate  Command  for  ACM  Entry  The  yaw  rate  command  is  simulated 
with  a  magnitude  of  1  degree,  second  for  two  seconds.  The  response  plots  for  all 
25  cases  can  be  found  in  Figure  G.38  through  Figure  G.45  of  appendix  G.  The 
surfaces  return  to  zero  after  the  command,  except  for  cases  6  and  7  which  establish 
new  equilibrium  positions  for  all  the  surfaces  of  roughly  1  and  4  degrees  from  the 
initial  trim  position.  The  rates  are  minimal.  Uncommanded  states  are  negligible 
and  return  to  zero.  The  commanded  state  is  identical  for  all  25  cases  but  the 
degree  of  sideslip  before  returning  to  zero  varies  by  roughly  0.5  degrees,  except 
again  for  cases  6  and  7.  These  two  cases  establish  a  constant  sideslip  of  about  2 
degrees  which  is  due  to  the  new  trim  position  of  the  surfaces. 

The  magnitude  of  the  second  yaw  rate  command  is  5  degrees /seconds.  The 
rudder  saturates  for  case  7  and  nearly  saturates  for  case  6.  All  rates  are  driven  to 
zero  except  case  7  which  starts  a  divergence  at  12  seconds  into  the  simulation.  The 
uncommanded  states  display  an  increased  transient  response  than  seen  previously 
for  the  small  magnitude  input.  Eventually,  the  uncommanded  states  then  return 
to  zero  for  almost  all  the  cases.  The  commanded  yaw  rate  is  identical  for  all 
simulations.  The  response  plots  for  all  25  cases  can  be  found  in  Figure  G.46 
through  Figure  G.53  of  appendix  G. 

Pitch  Rate  Command  for  TFTA  The  entire  set  of  single  input  commands 
was  to  be  simulated  for  the  TFTA  condition.  However,  upon  examination  of  the 
nominal  CRCA  response  to  a  2  degree /second  command,  further  simulations  are 
not  performed.  The  control  surfaces  displayed  a  high  frequency  flutter  with  a  small 
magnitude.  The  pitch  rate  response  in  figure  5.4  exhibits  a  high  frequency  mode 
overlaid  on  the  nominal  response  which  severely  distorts  the  response.  However,  the 
pitch  angle  response  is  quite  similar  to  that  achieved  in  the  ACM  entry  simulation. 
The  uncommanded  roll  and  yaw  rate  are  extremely  small  but  appear  to  diverge 
after  18  seconds  of  simulation. 
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Figure  5.4.  Pitch  Rate  Response  of  Nominal  CRCA  for  TFTA  Condition 

Conclusion  of  Single  Command  Input  Simulations  The  objective  of  the  single 
input  commands  is  to  demonstrate  decoupling  and  transient  response.  Overall, 
the  degree  of  cross-coupling  is  shown  to  be  quite  small  for  all  aircraft  within  the 
reduced  design  set  for  ACM  entry.  In  addition,  no  overshoot  or  oscillation  is  seen 
in  a  commanded  rate  state  or  associated  angle.  The  longitudinal  mode  is  shown 
to  be  the  most  robust  to  surface  failures  which  included  100/effectiveness  of  the 
left  canard.  The  pitch  rate  state  is  seen  to  have  the  fastest  rise  time  while  the 
two  lateral  rate  states  are  slightly  slower.  Minimal  lateral  acceleration  is  seen  for 
a  longitudinal  command.  However,  lateral  commands  develop  a  sizeable  amount 
of  longitudinal  acceleration.  The  negative  longitudinal  acceleration  for  the  pulse 
roll  command  can  be  partially  attributed  to  the  lack  of  thrust  modeling  in  the 
equations  of  motion.  With  the  aircraft  in  a  banked  orientation  and  no  thrust, 
generation  of  lift  is  reduced,  vertical  velocity  increases,  and  hence  the  negative 
acceleration. 

During  the  roll  simulations,  the  control  surfaces  did  not  return  to  their  equi- 


librium  positions  upon  completion  of  the  pulse  command.  While  the  rates  went 
to  zero,  the  surface  position  deflections  continued.  As  a  result,  the  uncommanded 
states  such  as  sideslip  are  seen  to  slowly  increase.  The  drift  after  the  command  is 
complete  is  attributed  to  the  use  of  rates  as  the  controlled  states.  In  all  cases,  the 
rates  are  driven  to  zero  after  the  command.  However,  once  the  rates  are  driven 
to  zero,  the  control  system  can  no  longer  regulate  the  associated  angles.  As  a 
consequence,  the  control  surfaces  begin  to  drift  in  response  to  the  dynamics  of  a 
nontrimmed  orientation. 

Simulations  for  the  TFTA  condition  are  not  performed.  The  constant  gain 
FCS  which  is  designed  for  aerodynamic  conditions  of  0.9  Mach  and  30000  feet  is 
unstable  for  aerodynamic  conditions  of  0.9  Mach  and  200  ft. 

Coordinated  Turn 

A  coordinated  turn  for  the  ACM  manuever  is  simulated  by  a  pulse  roll  com¬ 
mand  of  74.5  degrees /second  for  0.5  seconds,  a  step  pitch  rate  command  of  1.46 
degrees /second,  and  a  step  yaw  rate  command  of  1.46  degrees /second.  The  re¬ 
sponse  plots  for  all  25  cases  can  be  found  in  Figure  G.54  through  Figure  G.61  of 
appendix  G.  All  the  commands  are  ramped  as  in  the  previous  simulations.  For 
most  of  the  cases,  the  control  surfaces  move  rapidly  but  well  within  tolerances 
during  the  duration  of  the  pulse  roll  command.  The  surfaces  start  to  drift  after 
the  pulse  roll  command  as  observed  in  the  earlier  single  command  simulations  for 
roll  or  yaw  rates.  Cases  6,7,  and  to  a  smaller  extent  case  25  display  the  most  drift. 
All  control  surface  rates  are  driven  uniformly  to  zero  except  the  rudder:  several 
damage  cases  exhibited  roughly  a  0.5  degrees  second  rate  in  the  steady  state.  The 
three  commanded  rates  are  uniform  for  all  members  in  the  design  set.  The  bank 
angle  and  pitch  angle  are  indistinguishable  for  all  aircraft  within  the  reduced  de¬ 
sign  set.  However,  a  greater  range  of  sideslip  is  seen  with  an  average  of  about  2 
degrees  for  all  cases  except  cases  6  and  7. 
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Fourth  Order  Actuators 


The  first  order  actuators  are  replaced  by  the  fourth  order  actuators  given  in 
equation  2.9.  Without  any  changes  in  the  digital  controllers,  the  closed  loop  system 
is  unstable.  A  digital  filter  is  developed  which  is  cascaded  with  each  controller.  The 
numerator  of  the  filter  is  the  z  plane  mapping  of  the  three  new  poles  introduced  by 
the  fourth  order  actuators.  The  denominator  consists  of  three  nondominant  poles. 
The  resulting  filter  is  given  by 


F(z) 


0.39(2  -  0.42667  -  0.29457)(2  -  0.23841) 
(2  -  0.5)3 


(5.1) 


The  purpose  of  the  digital  filter  is  to  restore  the  original  loop  transmission  char¬ 
acteristics  by  cancelling  the  three  new  s  plane  poles  of  the  fourth  order  actuator 
with  the  three  2  plane  zeros.  The  three  nondominant  poles  do  not  add  any  phase 
or  magnitude  information  in  the  frequency  range  of  interest.  If  the  original  loop 
transmission  characteristics  are  restored,  system  performance  similar  to  the  origi¬ 
nal  system  with  first  order  actuators  is  expected. 

The  first  simulation  is  performed  with  a  six  degrees/second  pulse  command 
for  three  seconds.  The  magnitudes  of  the  surface  deflections  and  rates  for  the  pitch 
rate  simulation  are  roughly  proportional  to  that  obtained  with  2  degree/second 
pitch  command  with  first  order  actuators.  However,  the  three  damage  cases  sim¬ 
ulating  25%,  50%,  and  100%  left  canard  loss  exhibit  a  damped  oscillatory  motion 
of  the  control  surfaces  not  seen  in  the  earlier  simulations.  In  addition,  the  states 
exhibit  a  larger  variation  in  responses  between  the  bracketing  cases.  Where  the 
simulations  with  first  order  actuators  display  identical  pitch  rate  state  responses  for 
all  25  cases,  the  fourth  order  actuators  simulation  produce  state  responses  which 
do  not  track  the  pulse  command  well  in  several  cases.  In  addition,  where  the  un¬ 
modified  simulations  produced  a  2  degrees /second  response  for  a  2  degrees /second 
command,  the  modified  actuat  ors  produced  roughly  8  degrees  /second  for  a  6  de¬ 
grees/second  command.  Attempts  to  further  limit  the  error  by  varying  the  gain 


on  the  digital  filter  resulted  in  an  unstable  system.  Uneommanded  rate  states  and 
associated  angles  are  maintained  near  zero  but  their  transient  behavior  is  much 
more  oscillatory  than  seen  in  the  earlier  unmodified  simulation.  The  response  plots 
for  all  25  cases  can  be  found  in  Figure  G.62  through  Figure  G.69  of  appendix  G. 

A  second  simulation  is  performed  with  a  90  degrees /second  pulse  pitch  com¬ 
mand  for  0.5  seconds.  The  surface  deflections  show  a  damped  oscillatory  behavior 
which  is  similar  to  the  responses  achieved  with  a  first  order  actuator.  Cases  6,7, 
and  25  again  display  large  surface  movements  but  with  an  oscillatory  mode  over¬ 
laid  on  the  surface  excursion.  All  surface  rates  are  driven  to  zero  after  a  series 
of  medium  damped  oscillation  except  for  two  rudder  cases.  The  rudder  for  cases 
6  and  7  saturates  much  quicker  than  in  the  first  order  actuator  simulations.  The 
roll  rate  responses  display  some  oscillation  after  the  command  but  are  maintained 
at  zero.  The  response  plots  for  all  25  cases  can  be  found  in  Figure  G.70  through 
Figure  G.77  of  appendix  G. 

Overall,  better  system  response  is  obtained  in  the  longitudinal  mode  than 
in  the  lateral  mode  with  the  fourth  order  actuators.  In  neither  case  can  the  loop 
transmission  be  said  to  be  completely  restored.  Oscillatory  behavior  is  seen  in 
the  surface  deflections,  rates,  and  the  states  which  is  not  present  in  the  earlier 
simulations.  The  lateral  mode  oscillations  are  of  a  higher  frequency  than  the 
longitudinal  mode  oscillations.  The  oscillation  is  seen  in  almost  all  the  lateral 
states.  In  most  cases,  the  commanded  state  responses  of  pitch  rate  and  roll  rate 
are  seen  to  be  quite  similar  to  their  responses  obtained  with  first  order  actuators. 

Examination  of  Weighting  Matrix  Output 

A  final  simulation  is  performed  to  examine  the  individual  components  which 
form  the  pitch  rate  output.  Two  cases  are  examined:  the  nominal  case  and  and 
the  100%  canard  failure  case  using  both  first  order  and  fourth  order  actuators.  Ex¬ 
amination  of  the  transfer  functions  found  in  appendix  C  reveals  that  the  pitch  rate 
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state  can  be  viewed  as  the  composite  of  the  outputs  of  the  eight  open  loop  trans¬ 
fer  functions  (with  the  appropriate  feedback).  Each  transfer  function  represents 
a  portion  of  the  total  system  response.  The  weighting  matrix  is  used  to  convert 
a  single  input  into  eight  inputs  with  different  relative  authorities.  The  weighting 
matrix  outputs,  via  the  actuator  models,  are  then  the  inputs  to  the  eight  transfer 
functions  which  make  up  the  pitch  rate  output.  The  objective  of  the  weighting 
matrix  is  to  maximize  the  degree  to  which  the  surfaces,  and  therefore  the  control 
inputs  to  the  plant,  work  together.  If  the  control  surfaces  work  against  each  other, 
then  greater  surface  deflections  will  be  necessary  to  implement  a  particular  com¬ 
mand.  A  study  of  the  individual  component  outputs  is  an  indication  of  the  degree 
the  control  surfaces  begin  to  fight  each  other  under  simulated  failure  conditions. 

Examination  of  the  outputs  of  the  eight  pitch  rate  transfer  functions  by  them¬ 
selves  is  not  possible  since  each  has  an  unstable  pole.  Therefore,  an  alternative 
method  is  devised.  The  time  histories  of  the  outputs  of  the  eight  actuators  are  ob¬ 
tained  for  a  2  degree/second  pitch  rate  command.  The  actuator  outputs  are  taken 
because  the  weighting  matrix  outputs  undergo  an  integration  by  the  actuators  be¬ 
fore  going  to  the  plant.  Each  time  history  is  multiplied  by  the  its  corresponding 
value  in  the  third  row  of  the  input  B  matrix  from  chapter  two;  first  time  history 
multiplied  by  the  B(3,l)  term,  second  by  B(3,2),  and  so  on.  The  time  histones  are 
obtained  for  the  healthy  aircraft  and  the  100%  canard  loss  as  are  the  corresponding 
input  matrices  for  the  two  cases.  The  results  are  plotted  in  terms  of  symmetrical 
control  surfaces  and  can  be  found  in  Figures  5.5  through  5.7. 

The  results  of  the  simulation  are  not  unexpected  in  the  healthy  aircraft. 
The  symmetrical  surfaces  exert  equal  amounts  of  effort  and  work  together.  The 
failure  case  is  seen  to  display  unequal  efforts  in  terms  of  the  magnitude  of  the 
control  surface  motions.  For  the  failure  case,  the  canards  displayed  slightly  different 
responses  but  all  six  trailing  flaps  displayed  different  levels  of  response.  However, 
no  surfaces  are  seen  to  work  against  each  other  by  moving  in  opposite  directions. 
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The  first  order  actuators  are  replaced  by  fourth  order  actuators.  Examination 
of  the  nominal  aircraft  again  shows  the  surfaces  working  together  and  with  equal 
effort.  The  failure  case  clear  displays  conflicting  control  inputs  from  symmetrical 
surfaces.  Opposing  inputs  from  symmetrical  control  surfaces  will  generate  a  rolling 
and/or  yawing  moment  which  might  account  for  some  of  the  inertial  cross-coupling 
observed  for  failure  cases.  As  an  aside,  the  conflicting  surface  movement  is  not  seen 
for  all  surface  failure  cases;  a  simulation  run  with  a  100%  outer  flap  loss  did  not 
manifest  any  conflicting  surface  inputs. 
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Figure  5.5.  Control  Inputs  for  Unimpaired  CRCA  With  First  Order  Actuators 
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(  Figure  5.6.  Control  Inputs  for  CRCA  With  100%  Left  Canard  Loss  and  First  Order 
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Conclusion 

The  closed  loop  aircraft  is  simulated  with  a  variety  of  inputs.  Cross-coupling 
between  lateral  and  longitudinal  states  is  seen  to  be  held  to  a  minimum.  Cou¬ 
pling  between  the  lateral  states  for  a  lateral  command  is  somewhat  larger  in  some 
failure  cases.  The  objective  of  reducing  overshoot  in  pitch  rate  is  obtained  with 
first  order  actuators.  The  method  of  substituting  of  fourth  order  actuators  into  a 
system  designed  for  first  order  actuators  does  not  prove  to  be  wholly  satisfactory. 
Restoration  of  the  loop  transmission  characteristics  is  not  complete.  Given  suffi¬ 
cient  time,  improvement  of  the  loop  characterstics  may  be  possible  with  improved 

pole  locations  in  the  digital  filter.  Examination  of  the  weighting  matrix  outputs 

/ 

reveals  conflicting  control  inputs  for  at  least  one  failure  case. 
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VI.  Conclusions  and  Recommendations 


Conclusions 

The  application  of  QFT  to  a  sampled  data  system  is  shown  to  be  success¬ 
ful.  The  aircraft  FCS  continues  meeting  performance  objectives  in  spite  of  control 
surface  damage  in  the  simulations.  Pitch  overshoot  is  nonexistent  for  a  pitch  rate 
command  and  minimal  inertial  cross-  coupling  occurs  for  the  aircraft  within  the 
reduced  design  set.  Damped  oscillations  are  not  observed  in  any  of  the  controlled 
variables  or  associated  angles.  Once  the  weighting  matrix  is  derived,  the  design 
process  is  highly  transparent  with  the  trade-offs  clearly  detailed.  The  surface  de¬ 
flections  are  seen  in  the  simulations  to  be  minimal.  The  design  is  tolerant  to  control 
surface  failures.  In  conclusion,  the  use  of  QFT  to  design  a  FCS  for  an  unstable, 
nonminimum  phase  aircraft  is  effective  for  the  conditions  upon  which  the  design  is 
based. 

However,  extensions  to  conditions  other  than  the  design  conditions  such  as 
the  TFTA  aerodynamic  condition  or  the  use  of  fourth  order  actuators  is  not  as 
successful.  The  use  of  the  fourth  order  actuators  does  not  have  a  large  impact 
on  the  response  of  longitudinal  mode  states  to  a  pitch  rate  command.  However, 
the  lateral  mode  responses  to  a  roll  rate  command  are  divergent  for  several  failure 
cases.  The  design  is  unacceptable  for  direct  application  to  the  TFTA  condition. 
The  lack  of  success  in  extending  the  design  beyond  the  original  conditions  only 
emphasizes  the  need  to  incorporate  such  information  into  the  design  process  from 
the  beginning. 

The  following  discussions  highlight  some  of  the  issues  raised  during  the  design 
process. 

Weighting  Matrix  Different  values  other  than  zero  are  originally  used  in  the 
off-diagonal  terms  of  the  left  hand  side  of  equation  (2.12).  The  objective  is  to 


mirror  physical  reality  in  that  complete  decoupling  between  states  is  not  possible. 
However,  this  approach  is  discarded  because  the  terms  of  the  resulting  weighting 
matrix  no  longer  made  physical  sense.  The  magnitudes  of  terms  in  the  weighting 
matrix  for  symmetrical  surfaces  such  as  left  and  right  canard  are  found  to  be 
unequal.  Different  weightings  for  similar  surfaces  would  lead  to  an  untrimmed 
aircraft.  As  the  off-diagonal  terms  of  the  left  hand  side  of  equation  (2.12)  decrease 
in  magnitude,  the  difference  between  symmetrical  terms  of  the  weighting  matrix 
become  smaller.  Zeros  are  therefore  used  in  the  off-diagonal  terms  and  the  resulting 
weighting  matrix  is  symmetrical  with  respect  to  like  surfaces. 

While  the  weighting  matrix  is  developed  in  the  continuous  plane,  the  mini¬ 
mum  phase  properties  of  the  effective  plant  are  verified  in  the  discrete  wi  plane.  As 
an  intermediate  step  in  the  conversion  from  the  continuous  to  the  discrete  domain, 
the  effective  plant  is  modified  by  cascading  the  actuator  model  with  each  transfer 
function  in  the  plant  matrix.  When  the  determinant  of  the  effective  plant  is  ex¬ 
amined  in  the  s  plane  without  cascading  the  actuator  model,  additional  cases  of 
nonminimum  phase  plant  matrices  are  found.  The  correlation  between  the  weight¬ 
ing  matrix  and  when  best  to  determine  whether  the  plant  is  minimum  phase  or 
not  is  not  clear.  Further  study  is  needed  to  determine  if  the  test  for  an  adequate 
weighting  matrix  should  be  performed  in  the  continuous  or  discrete  domain  and 
with  or  without  the  actuators  cascaded  with  the  basic  plant. 

The  use  of  the  weighting  matrix,  while  not  truly  a  part  of  QFT,  is  imposed  by 
the  requirements  of  QFT.  The  method  used  to  develop  the  terms  of  the  weighting 
matrix  appears  satisfactory.  However,  there  is  no  indication  at  the  onset  of  the 
design  process  if  the  weighting  matrix  derived  in  this  thesis  is  the  optimal  weighting 
matrix  for  the  CRCA.  The  examination  of  the  component  responses  of  pitch  rate 
output  reveals  under  some  failure  conditions  the  surfaces  are  fighting  each  other. 
Obviously,  optimal  performance  under  failure  conditions  can  not  be  achieved  if  the 
control  inputs  are  not  compatible  with  each  other.  The  full  benefits  of  QFT  may 


>-2 


not  be  realized  on  nonsquare  plants  if  the  weighting  matrix  is  developed  without 
some  indications  on  how  to  minimize  the  amount  of  crossed  control  inputs. 

TFTA  The  control  system  designed  for  the  ACM  entry  condition  could  not 
be  extended  to  include  the  TFTA  condition  without  some  type  of  gain  scheduling. 
The  difference  between  the  two  conditions  is  the  relative  aerodynamic  pressure. 
The  design  philosophy  applied  in  this  thesis  is  to  maximize  the  loop  transmission. 
Therefore,  the  design  is  optimized  for  the  aerodynamic  conditions  represented  by 
ACM  entry.  A  possible  approach  to  extending  the  design  to  TFTA  would  be  to 
reduce  the  loop  transmission.  Examination  of  the  loop  transmission  plot  similar 
to  the  one  in  Figure  4.1  reveals  the  upper  and  lower  limits  on  the  allowable  gain 
adjustments  before  the  stability  boundaries  are  violated.  The  probable  tradeoff 
would  be  to  sacrifice  performance  at  the  ACM  entry  condition  while  increasing 
performance  in  the  TFTA  condition. 

Fourth  Order  Actuators  The  QFT  design  is  based  upon  the  cascade  com¬ 
bination  of  the  effective  plant  and  the  first  order  actuators.  After  the  z  plane 
controller  is  implemented,  a  simulation  is  performed  using  fourth  order  actuators 
without  changing  the  controller.  The  replacement  causes  the  system  to  become 
unstable.  The  two  actuator  models  have  very  similar  step  responses  in  the  time 
domain  which  would  initially  indicate  the  two  models  could  be  interchanged.  An 
explanation  for  the  unstable  response  lies  in  the  frequency  domain.  The  phase 
response  of  the  two  models  is  quite  different.  The  first  order  model  contributes  a 
maximum  of  -90  degrees  of  phase  angle  while  the  fourth  order  model  will  contribute 
up  to  -360  degrees.  The  difference  in  phase  angle  is  critical  when  one  considers  that 
the  loop  transmission  is  shaped  based  on  a  desired  loop  phase  response.  There¬ 
fore,  any  changes  in  the  phase  response  will  obviously  have  a  dramatic  effect  in  the 
regions  where  the  stability  boundaries  are  rigorously  satisfied. 

The  use  of  the  digital  filter  with  the  fourth  order  actuators  almost  restores 
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the  loop  transmission  characteristics  of  the  system  with  first  order  actuators.  The 
system  performance  for  a  pitch  rate  command  is  quite  similar  to  that  achieved  with 
the  first  order  actuators  except  for  the  increased  magnitude  of  the  commanded  state 
and  the  oscillatory  behavior  of  the  states  and  control  surfaces.  The  response  of  the 
aircraft  for  the  roll  rate  command  is  not  satisfactory.  Several  failure  cases  diverged 
with  fourth  order  actuators  that  are  stable  with  the  first  order  actuators.  One 
conclusion  drawn  is  the  use  of  first  order  actuators  obscures  information  which  is 
then  unaccounted  for  in  the  design.  The  poor  performance  of  the  roll  rate  command 
simulation  indicates  the  unaccounted  information  is  critical  and  needs  to  be  built 
into  the  design. 

The  gain  on  the  filter  is  originally  such  that  the  steady  state  error  of  the 
filter  is  zero.  However,  the  closed  loop  system  is  unstable  so  the  gain  is  varied 
until  the  desired  performance  is  obtained.  The  best  performance  is  obtained  with 
the  gain  equal  to  0.39;  a  lower  value  caused  the  magnitude  of  the  response  to 
increase  dramatically  while  a  larger  gain  value  induced  oscillation  and  eventually 
instability. 

The  approach  of  designing  with  first  order  actuators  and  simulating  with 
fourth  order  actuators  is  imperfect  at  this  stage.  Insufficient  time  remained  to 
fully  explore  alternative  pole  postions  or  gains.  The  use  of  the  filter  in  the  discrete 
domain  to  cancel  the  poles  of  the  actuators  in  the  continuous  domain  has  the 
unintended  effect  of  amplifying  the  command  response.  In  addition,  the  oscillatory 
behavior  due  to  the  set  of  complex  poles  in  the  fourth  order  actuators  is  not 
adequately  controlled  especially  in  the  lateral  mode.  On  the  positive  side,  the 
technique  does  allow  the  designer  to  study  the  overall  system  behavior  with  higher 
order  actuators.  In  the  case  of  this  thesis,  the  increased  oscillatory  behavior  seen 
with  the  higher  order  actuators  indicates  a  shortcoming  of  designing  with  simplified 
actuators.  Such  results  only  further  stress  the  need  to  develop  techniques  which 
incorporate  the  higher  order  actuators  into  the  plant  model  before  the  loop  shaping 


is  performed. 
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Digital  Aspects  The  use  of  the  u>/  plane  to  carry  out  the  design  appears  to 
be  an  excellent  alternative  to  designing  in  the  z  plane.  The  use  of  the  Hoffman 
algorithm  allows  the  equations  of  motion  to  be  transformed  from  the  s  plane  to  the 
wi  plane  without  manipulating  the  intermediate  results  in  the  z  plane.  Once  in  the 
wi.  the  loop  transmission  function  is  shaped  based  on  Professor  Horowitz's  repeated 
exhortations  to  ‘'pack  in  as  much  loop  transmission  as  possible’’.  The  reason 
behind  this  approach  is  based  on  the  inherent  limitation  of  sampling.  However, 
this  approach  results  in  a  high  open  loop  bandwidth  in  the  vicinity  of  100  rad /sec 
in  the  w>  plane.  Time  did  not  permit  a  study  of  the  correlation  between  a  high 
open  loop  bandwidth  in  the  wt  plane  and  the  actual  bandwidth  in  the  z  plane. 
With  the  increasing  use  of  digital  flight  control  systems,  an  investigation  into  the 
frequency  relationships  of  the  s.  z,  and  wt  is  warranted. 

Recommendations 

1.  A  digital  QFT  control  system  with  a  sample  rate  of  80  Hz  for  the  CRCA 
should  be  performed  with  tracking  and/  or  disturbance  bounds.  The  stability 
bounds  in  this  thesis  are  easily  satisfied  and  therefore  the  further  constraint  of 
tracking  and/or  disturbance  bounds  should  be  within  reason.  Once  the  design 
is  complete,  further  thesis  efforts  should  be  directed  at  fully  characterizing 
the  design.  Susceptibility  to  noise  inputs  from  the  sensors  or  in  terms  of 
wind  gusts  and  man-in-the-loop  simulations  are  examples  of  areas  of  possible 
exploration.  The  objective  is  to  demonstrate  the  qualitative  aspects  of  a 
QFT  based  flight  control  system. 

2.  The  use  of  the  CRCA  equations  of  motion  provide  an  excellent  opportu¬ 
nity  for  further  thesis  efforts  regardless  of  technique  applied.  The  aircraft  is 
representative  of  an  advanced  airframe  with  its  attendant  problems  of  insta- 
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bilities  and  nonminimum  phase  characteristics.  If  flight  control  research  at 
AFIT  is  to  continue  to  be  relevant  to  outside  agencies,  research  efforts  should 
concentrate  on  the  types  of  problems  presented  by  the  CRCA. 

3.  The  software  package  MatrixX  is  an  excellent  method  to  simulate  dynamic 
systems.  However,  there  are  limits  to  its  usefulness  as  a  design  package  due 
to  numerical  problems.  Limitations  exist  within  MatrixX  on  the  number  of 
variables  which  can  be  used,  the  accuracy  available  in  the  manipulations  of 
complex  numbers,  the  maximum  order  of  polynomials  which  can  be  handled, 
and  reliability  with  which  phase  angles  are  calculated  for  unstable  or  non¬ 
minimum  phase  systems.  The  numerical  limitations  of  MatrixX  need  to  be 
understood  and  disseminated  to  future  AFIT  students. 

4.  The  development  of  the  terms  within  the  weighting  matrix  need  to  be  corre¬ 
lated  with  some  measure  of  optimality.  Until  the  weighting  matrix  is  shown 
to  be  optimal,  the  full  benefits  of  QFT  on  nonsquare  plants  may  not  be  re¬ 
alized.  For  instance,  an  optimal  weighting  matrix  may  eliminate  the  small 
amount  of  sideslip  observed  for  the  two  failure  cases  which  differed  markedly 
from  the  other  23  cases  in  the  pitch  rate  simulations.  Examination  of  the 
individual  component  responses  for  a  particular  output  reveals  at  least  part 
of  the  unwanted  inertial  cross-coupling  may  be  due  to  control  surfaces  work¬ 
ing  against  each  other.  The  weighting  matrix  for  this  thesis  is  developed 
based  on  the  healthy  aircraft.  An  alternative  approach  might  be  to  develop 
the  weighting  matrix  based  on  the  worst  case  aircraft.  Another  possibility 
is  to  average  the  terms  of  several  weighting  matrices  developed  for  several 
different  types  of  failures  to  produce  a  composite  weighting  matrix. 

5.  The  states  used  for  feedback  in  this  thesis  are  pitch  rate,  roll  rate,  and  yaw 
rate.  Once  the  rates  are  driven  to  zero,  control  over  the  associated  angles 
is  lost.  The  situation  is  especially  exaggerated  in  the  rolling  mode  since  an 
intrinsic  re^+oring  moment  does  not  exist.  The  conclusion  drawn  is  the  use 
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of  only  rates  is  not  a  satisfactory  set  of  feedback  states.  A  better  set  might 
be  pitch  rate,  roll  rate,  and  sideslip. 

6.  The  relationship  between  higher  order  actuators  and  nonminimum  phase 
zeros  in  the  discrete  domain  should  be  investigated.  Techniques  need  to 
be  developed  that  allow  the  incorporation  of  higher  order  actuators  into  the 
basic  plant  before  loop  shaping  is  accomplished.  As  part  of  this  investigation, 
the  restoration  of  the  loop  transmission  characteristics  through  the  use  of  the 
digital  filter  should  be  further  explored. 
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Appendix  A.  Tracking  Response  Models 


The  tracking  response  models  are  taken  from  previous  theses  which  dealt 
with  a  QFT  design  for  the  AFTI  Fl6.  The  response  models  are  given  first  for  the 
continuous  case  and  then  for  the  wt  case.  The  wt  versions  are  obtained  by  using 
the  Hoffman  algorithm  to  transform  from  the  5  plane  through  the  z  plane  to  the 
wt  plane. The  sample  rate  is  80  Hz. 


Tu(s)  = 


Roll  Rate  Model 

Upper  Bound 
2.136(.s  -  1.3) 


s2  +  2.667s  -  2.776 


Tu{rvt)  = 


—0.013 3(w2t  -  158.7tcf  -  208) 
ti'2r  -  2.667 wt  -  2.776) 

Lower  Bound 
6.937 


r,(u*o  = 


T‘(3'  (s  —  0.75)(s  +  2.5)(s  -r  3.7) 

5.642  •  10_7(tr  -  160)  (tc  -  280.7)  (tc  -  273.7) 
(u-  -  3.699)  (tc  -  2.5)(tc  -  0.75) 

Pitch  Rate  Model 

Upner  Bound 
4.167(s  +  2.667) 


T»(s)  = 


s2  +  5.334s  +  11.113 


= 


-2.604  •  10_2(tc  -  160) (tw  •+■  2.667) 
(tc  +  2.667  ±  j'2.0) 

Lower  Bound 


TM  = 


63.75 


Tu(wt)  = 


(s  +  1.5)(s  +  4.25)(s  +  10) 

5.173  •  10~6(tc  -  160) (tc  -  269.5) (tc  -  285.3) 
(tc  +  9.987)  (u-  +  4.249)  (tc  +  1.5) 


1 


■ 


1 


| 
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Yaw  Rate  Model 

Upper  Model 


4.167(.s  -f-  2.667) 

52  -  5.334s  -r  11.113 


Tu(wi) 


-2.604  •  10- 2 (w  -  160) (u,'  +  2.667) 
(u/  +  2.667  ±  j2.0) 

Lower  Model 


Tl ^  (s  -r  1.5)(s  4.25)(s  -  10) 

5.173  •  10-6(u'  -  160)(tr  -  269.5)  (tr  -  285.3) 
(tr  +  9.987)(u>  -r  4.249)(tr  -  1.5) 
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Appendix  B.  Control  Surface  Failure  List 

Case  Failure  Description  for  ACM  Entry 

1  unimpaired  aircraft 

2  25%  rudder  loss 

3  50%  rudder  loss 

4  100%  rudder  loss 

5  25%  outer  flap  loss 

6  50%  outer  flap  loss 

7  100%  outer  flap  loss 

8  25%  inner  flap  loss 

9  50%  inner  flap  loss 

10  100%  inner  flap  loss 

11  25%  outer  flap  and  25%  elevator  loss 

12  50%  outer  flap  and  50%  elevator  loss 

13  100%  outer  flap  and  100%  elevator  loss 

14  25%  elevator  loss 

15  50%  elevator  loss 

16  100%  elevator  loss 

17  elevator  fail  to  trail 

18  canard  locked  at  5  degrees 

19  rudder  fail  to  trail 

20  canard  fail  to  trail 

21  outer  flap  fail  to  trail 

22  inner  flap  fail  to  trail 

23  canard  and  rudder  fail  to  trail 

24  canard  and  outer  flap  fail  to  trail 

25  outer  flap  and  elevator  fail  to  trail 

26  25%  canard  loss 

27  50%  canard  loss 

28  100%  canard  loss 

The  reduced  design  set  for  ACM  entry  consists  of  cases  1  to  19,21,23,25  to  28. 
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Case 


Failure  Case  List  for  TFTA 


1 

2 

8 

4 

5 

6 

7 

8 

9 

10 
11 
12 


13 

14 

15 

16 

17 

18 


unimpaired  aircraft 
canard  and  outer  flap  fail  to  trail 
rudder  fail  to  trail 
outer  flap  fail  to  trail 
canard  fail  to  trail 
25*3?  outer  flap  loss 
50%  outer  flap  loss 
100%  outer  flap  loss 
25%  outer  flap  and  25%  elevator  loss 
50%  outer  flap  and  50%  elevator  loss 
100%  outer  flap  and  100%  elevator  loss 
outer  flap  and  elevator  fail  to  trail 
25%  rudder  loss 
50%  rudder  loss 
100%  rudder  loss 
25%  canard  loss 
50%  canard  loss 
100%  canard  loss 
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Appendix  C.  Plant  Open  Loop  Transfer  Functions  for  Nominal 

ACM  Case  in  the  s  Plane 

Actuator  dynamics  are  not  included  in  the  following  transfer  functions. 

Pitch  Rate  Transfer  Functions 

0.102s(s  -  1.035)(s  -  0.0096) 

Pl1  ”  (s  -  1.694) (s  -  3.37)  (s  -  0.00402  =  >0.0356) 

0_  0.102s(fi  -  1.035)(s  -  0.0096) 

P  ~  (s  -  1.694)  (s  ~  3.37) (s  -  0.00402  -  j0.0356) 

-0.028 s(s  -  1.296)(s  +  0.0102) 

Pl  “  (s  -  1.694)(s  +  3.37)  (5  +  0.00402  ±  j0.0356) 

-0.022s{s  -  1.253)(s  +  0.0101) 

P  ~  (s-  1.694) (s  4-  3.37) (s  +  0.00402  ±  ;0.0356) 

-0.028s(s  -  1.296)  (s  -  0.0102) 

P  “  (s  -  1.694) (s  -  3.37) (s  -  0.00402  -  j0.0356) 

-0.022s(s  -  1.253)(s  -  0.0101) 

P  “  (s  -  1.694)(s  *  3.37)(s  -r  0.00402  ±  ;0.0356) 

-0.02s(s  +  1.3l2)(s  +  0.0102) 

P  “  (s  -  1.694) +  3.37) (s  -f  0.00402  ±  j0.0356) 

Q_  -0.02s(s  1.312)(s  +  0.0102) 

P  ~  (5  -  1.694) (s  +  3  37)  (s  -  0.00402  ±  j0.0356) 


pl9  =  0 
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Roll  Rate  Transfer  Functions 


o  _  0.076(5  -0.00126)  (5-  0.4247  x  j4.505) 

P  1  "  {s  +  2.017)(s  ■+■  0.0562) (5  2.492  x  >1.834) 


-0.076(s  -  0.00126)(5  -  0.4247  1  ;4.505) 
P22  "  (s  ~  2.017)(«  -  0.0562)  (5  -  2.492  ±  >1.834) 


0.222(s  -  0.00124) (5  x  0.2695  *  >1.675) 
p23  '  (s  -  2.017)(s  -  0.0562) (5  ■+  2.492  -  >1.834) 


0.201(5  -  0.00124) (5  ■+•  0.2686  ±  >1.645) 

P  4  "  (s  +  2.017)(s  x  0.0562) (5  +  2.492  x  >1.834) 


_  -0.222(5  -  0.00124)(5  -  0.2695  ±  >1.675) 
P  ~  (5  +  2.017) (5  x  0.0562) (5  -  2.492  x  >1.834) 


_  — 0.20l(s  -  0.00124)(s  -  0.2686  x  >1.645) 
P  ~  (5  -  2.017)(5  +  0.0562) (5  -  2.492  x  >1.834) 


_  0.111(5  -  0.00124)(s  x  0.2706  x  >1.712) 

P  ~  (5  +  2.01 7) (5  -1-  0.0562) (5  -  2.492  x  >1.834) 


oo  __  -0.111(5  -  0.00124)(5  +  0.2706  ±  j  1.712) 
P  ~  (s  +  2.017)(s  +  0.0562) (s  +  2.492  ±  >1.834) 


0.115(s  +  3.298)(s  -  3.239) (5  -  0.00125) 

P  ~  (5  x  2.017)(s  1-  0.0562)(5  +  2.492  ±  >1.834) 
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Appendix  D.  Effective  pu(wt)  Open  Loop  Transfer  Functions  for 

Nominal  ACM  Entry  Case 


r 

f  ■ 

; 

Continuous  actuator  dynamics  of  (57.3) (7^5)  are  cascaded  with  each  plant  transfer 
r  functions  after  incorporation  of  the  weighting  matrix.  The  resulting  effective  plant 

transfer  functions  are  transformed  to  the  wt  plane.  The  sample  rate  is  80  Hz. 

_  1.4916  •  10-3u>/(u>/  +  3723) (wt  -  0.00969) (wt  -r  1.07) (ie/  -  160) 

Pn'U’  ~  {wi  -  19. 9){u>f  -  3.37)(u>f  -  1.6494)(u»  -  0.00402  -  >00356^ 


pl2{wi)  =  0 


pi3{wt)  =  0 


■ 


Pi  i(wf)  =  0 


p22(wr) 


-1.1595  ■  10“3(^  +  3491) (u;/  -  0.00125)(u>/  -  160)(u7  -  0.2763  ±  >1.890) 
(tut  +  19.9) (u;  +  2.017)(utf  +  0.0562)(uU  -  0.2492  it  >1.834) 


P23(^'0  = 


-2.1832  •  10~3(u;f  -  3539)(u>f  -  0.00125)(u;f  -  160)(u>f  +  0.4266  ±  3.989) 
(wt  +  19.9 )(wt  -r  2.017 )(wf  -h  0.0562)  (wt  7  0.2492  =  jl.834) 


p3j(u;r)  =  0 


i- 


P32  (utf)  = 


-4.5543  •  10~s(u>f  +  3657) (wtl.  150) {wf  -  160 )(u;f  +  0.2002  ±  >1.645) 
(wt  +  19.9)  (tuf  +  2.017)(u;t  +  0.0562)(u;f  +  0.2492  ±  >1.834) 


P33(W) 


-1.0022  •  10-3(u;f  +  3750)(wt  +  1.7l0)(u>/  -  160)(tef  +  0.1785  ±  >0.8124) 
{wt  +  19.9)(wt  +  2.017)(u;t  +  0.0562)(u;/  +  0.2492  ±  >1.834) 
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Appendix  E.  Effective  qu(w/)  Open  Loop  Transfer  Functions  for 

Nominal  A  CM  Entry  Case 


Continuous  actuator  dynamics  of  (57. 3)  (^_)  are  cascaded  with  the  plant  transfer 
functions  in  the  s  plane  before  transformation  to  the  wt  plane. 

1.4916  •  I0~3wt(wt  -  3723) (wt  +  0.00969) (wt  -  1.070)(u>f  -  160) 

9,1  ~  (wt  -  19.9 )(wt  -  3.369)  (tr/  -  1.6494) (trf  -  0.004018  =  >0.03561) 


qn(w>)  -  oc 


9n(u>»)  =  oc 


q2i(wf)  =  oc 


.  1. 0603  •  10-3u7(u,/ -  3489) (u;/ -  0.0628-  >0.0007) (wt  -  160) 

Q27\1L}t)  — - - - - — — - * 

(wt  19.9)(u,7  —  0.0562)(ir  +  1.713)(wf  +  0.1785  i  >0.8124) 

.  .  -4.867  •  10 ~4wt(wt  ~  3757) («/f  +  0.0628  ±  >0.0007) (wt  -  160) 

v  («;/  +  19.9)  (wt-  0.001255)  (wf  -  0.0562)  (wt  +  0.4266  ±>3.989) 


931  (iy/)  =  oc 


,  ,  -  2.3332 -i0-2u;/(u-f-  3489)(u’/-  0.0628  ij0.0007)(u^- 160) 

v  (wt  +  19.9)(u>/  —  1.15)(u;f  -  0.0562)(u^  0.2002  ±  >1.645) 

9.1645  •  10-4tyf(iyf  +  3757)(u;/  +  0.0628  ±  >0.0007)(u;f  -  160) 
9531^  )  -  ^  +  19  _  o.001245)(u;f  -  0.0562)(u»  +  0.2763  ±  >1.89) 
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Appendix  F.  Loop  Transmission  Characteristics 

Pitch  Rate  Loop 

The  following  table  displays  the  magnitude  and  phase  angle  obtained  for  the 
pitch  rate  loop  transmission  function. 


Frequency 

Magnitude  (dB) 

Phase  Angle  (degrees) 

0.001 

12.226 

-163.108 

0.002 

12.895 

-147.158 

0.004 

14.981 

-120.166 

0.008 

20.247 

-85.087 

0.016 

30.127 

-54.898 

0.032 

49.934 

-68.793 

0.064 

43.764 

-179.592 

0.128 

41.241 

-176.455 

0.256 

40.716 

-168.212 

0.512 

40.626 

-157.601 

1.024 

40.403 

-141.469 

2.048 

38.64 

-126.717 

4.096 

34.03 

-117.909 

8.192 

27.613 

-113.135 

16.384 

20.314 

-111.092 

32.768 

12.252 

-104.788 

65.536 

4.798 

-87.154 

131.07 

-0.559 

-62.128 

262.14 

-3.423 

-37.879 

524.29 

-4.486 

-21.294 

1048.6 

-4.797 

-12.999 

2097.2 

-4.939 

-11.110 

4194.3 

-5.426 

-14.599 

8388.6 

-6.900 

-18.423 

16777 

-8.601 

-15.386 

33554 

-9.475 

-9.3000 

67109 

-9.757 

-4.911 

134218 

-9.859 

0 

Roll  Rate  Loop 

The  following  table  displays  the  magnitude  and  phase  angle  obtained  for  the 
roll  rate  loop  transmission  function. 


Frequency 

Magnitude  (dB) 

Phase  Angle  (degrees) 

0.001 

42.805 

-177.314 

0.002 

42.807 

-174.629 

0.004 

42.817 

-169.270 

0.008 

42.855 

-158.640 

0.016 

43.005 

-138.040 

0.032 

43.561 

-101.371 

0.064 

45.265 

-50.153 

0.128 

48.715 

-5.006 

0.256 

53.202 

16.854 

0.512 

59.405 

17.995 

1.024 

62.992 

-118.365 

2.048 

46.577 

-138.503 

4.096 

35.387 

-138.976 

8.192 

25.529 

-131.623 

16.384 

16.724 

-122.700 

32.768 

8.127 

-110.832 

65.536 

0.471 

-89.186 

131.07 

-4.911 

-60.435 

262.14 

-7.559 

-31.793 

524.29 

-7.875 

-10.360 

1048.6 

-6.784 

-1.227 

2097.2 

-5.996 

-5.565 

4194.3 

-6.654 

-13.744 

8388.6 

-8.251 

-18.089 

16777 

-9.985 

-15.929 

33554 

-10.983 

-10.062 

67109 

-11.328 

-5.402 

134218 

-11.456 

0 
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The  following  table  displays  the  magnitude  and  phase  angle  obtained  for  the 
yaw  rate  loop  transmission  function. 


Frequency 

Magnitude  (dB) 

Phase  Angle  (degrees) 

0.001 

96.816 

-319.097 

0.002 

93.446 

-297.519 

0.004 

88.437 

-278.382 

0.008 

82.657 

-260.936 

0.016 

76.490 

-238.552 

0.032 

69.867 

-199.679 

0.064 

66.205 

-124.409 

0.128 

67.379 

-135.704 

0.256 

57.324 

-116.87 

0.512 

50.963 

-91.508 

1.024 

47.708 

-67.671 

2.048 

54.564 

-160.542 

4.096 

29.870 

-162.274 

8.192 

19.325 

-134.623 

16.384 

11.397 

-117.759 

32.768 

4.002 

-104.754 

65.536 

-2.550 

-90.188 

131.07 

-7.770 

-76.043 

262.14 

-11.801 

-67.715 

524.29 

-15.929 

-67.046 

1048.6 

-20.859 

-65.497 

2097.2 

-25.763 

-56.156 

4194.3 

-29.378 

-39.487 

8388.6 

-31.125 

-22.939 

16777 

-31.702 

-12.028 

33554 

-31.859 

-6.092 

67109 

-31.899 

-3.056 

134218 

-31.913 

0 

Appendix  G.  Simulation  Results 


A  CM  Entry  Condition 

Small  Magnitude  Pitch  Rate  Command  Response  The  commanded  pitch  rate 
has  a  magnitude  of  2  degrees /second  and  a  duration  of  three  seconds. 
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Figure  G.l.  Left  Canard  Position  and  Rate  Deflection  of  25  Cases  for  ACM  Entry 


G-l 


Figure  G.2.  Left  Outer  Flap  Position  and  Rate  Deflection  of  25  Cases  for  ACM 
Entry 


Figure  G.3.  Left  Inner  Flap  Position  and  Rate  Deflection  of  25  Cases  for  ACM 
Entry 
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Figure  G.5.  Rudder  Position  and  Rate  Deflection  of  25  Cases  for  ACM  Entry 
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Figure  G.6.  Pitch  Rate  Response  for  25  Cases  for  ACM  Entry 


Figure  G.7.  Pitch  Angle  (6)  for  25  Cases  for  ACM  Entry 
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Figure  G.8.  Roll  Rate  Response  for  25  Cases  for  ACM  Entry 


Figure  G.9.  Bank  Angle  (<t>)  for  25  Cases  for  ACM  Entry 
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Figure  G.ll.  Sideslip  (0)  for  25  Cases  for  ACM  Entry 
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Figure  G.13.  Lateral  Acceleration  for  25  Cases  for  ACM  Entry 
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Large  Magnitude  Pitch  Rate  Command  Response  The  commanded  pitch  rate 
■  has  a  magnitude  of  12  degrees/ second  and  a  duration  of  three  seconds. 


Figure  G.14.  Pitch  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.17.  Bank  Angle  (<i>)  for  25  Cases  for  ACM  Entry 
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Figure  G.18.  Yaw  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.20.  Normal  Acceleration  for  25  Cases  for  ACM  Entry 
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Small  Magnitude  Roll  Rate  Command  Response  The  commanded  roll  rate 


has  a  magnitude  of  30  degrees  second  and  a  duration  of  0.5  seconds. 


Figure  G.22.  Pitch  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.23.  Pitch  Angle  (0)  for  25  Cases  for  ACM  Entry 
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Figure  G.24.  Roll  Rate  Response  for  25  Cases  for  ACM  Entry 


Figure  G.27.  Sideslip  (/3)  for  25  Cases  for  ACM  Entry 
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Large  Magnitude  Roll  Rate  Command  Response  The  commanded  roll  rate 
has  a  magnitude  of  135  degrees  second  and  a  duration  of  0.5  seconds. 
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Figure  G.30.  Pitch  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.31.  Pitch  Angle  (0)  for  25  Cases  for  ACM  Entry 
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Figure  G.32.  Roll  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.33.  Bank  Angle  (0)  for  25  Cases  for  ACM  Entry 
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Figure  G.34.  Yaw  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.35.  Sideslip  (/?)  for  25  Cases  for  ACM  Entry 


Figure  G.36.  Normal  Acceleration  for  25  Cases  for  ACM  Entry 
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Small  Magnitude  Yaw  Rate  Command  Response  The  commanded  yaw  rate 
has  a  magnitude  of  1  degree  /second  and  a  duration  of  two  seconds. 


Figure  G.38.  Pitch  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.39,  Pitch  Angle  (0)  for  25  Cases  for  ACM  Entry 


Figure  G.40.  Roll  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.41.  Bank  Angle  (0)  for  25  Cases  for  ACM  Entry 


Figure  G.42.  Yaw  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.43.  Sideslip  (0)  for  25  Cases  for  ACM  Entry 


Figure  G.44.  Normal  Acceleration  for  25  Cases  for  ACM  Entry 
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Figure  G.45.  Lateral  Acceleration  for  25  Cases  for  ACM  Entry 
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Large  Magnitude  Yaw  Rate  Command  Response  The  commanded  yaw  rate 
has  a  magnitude  of  5  degrees  second  and  a  duration  of  two  seconds. 


Figure  G.46.  Pitch  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.47.  Pitch  Angle  (5)  for  25  Cases  for  ACM  Entry 


Figure  G.48.  Roll  Rate  Response  for  25  Cases  for  ACM  Entry 


Figure  G.49.  Bank  Angle  (<t>)  for  25  Cases  for  ACM  Entry 


Figure  G.50.  Yaw  Rate  Response  for  25  Cases  for  ACM  Entry 
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Figure  G.51.  Sideslip  {&)  for  25  Cases  for  ACM  Entry 
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Figure  G.52.  Normal  Acceleration  for  25  Cases  for  ACM  Entry 
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Coordinated  Turn 


The  commanded  pitch  rate  has  a  magnitude  of  2  degrees/second,  the  com¬ 
manded  yaw  rate  has  a  magnitude  of  2  degrees /second  and  the  commanded  roll 
rate  has  a  magnitude  of  2  degrees /second. 
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Figure  G.54.  Pitch  Rate  Response  for  25  Cases  for  Coordinated  Turn 
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Figure  G.55.  Pitch  Angle  (6)  for  25  Cases  for  Coordinated  Turn 
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Figure  G.56.  Roll  Rate  Response  for  25  Cases  for  Coordinated  Turn 


Figure  G.57.  Bank  Angle  ( <j> )  for  25  Cases  for  Coordinated  Turn 


Figure  G.58.  Yaw  Rate  Response  for  25  Cases  for  Coordinated  Turn 


Figure  G.59.  Sideslip  (3)  for  25  Cases  for  Coordinated  Turn 


Figure  G.60.  Normal  Acceleration  for  25  Cases  for  Coordinated  Turn 


Figure  G.61.  Lateral  Acceleration  for  25  Cases  for  Coordinated  Turn 
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Fourth  Degree  Actuators 

Pitch  Rate  Command  The  pitch  rate  command  has  a  magnitude  of  6  de¬ 
grees/second  and  a  duration  of  three  seconds. 
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Figure  G.62.  Pitch  Rate  Response  for  25  Cases  with  Fourth  Degree  Actuators 


G-38 


Figure  G.63.  Pitch  Angle  (0)  for  25  Cases  with  Fourth  Degree  Actuators 


Figure  G.64.  Roll  Rate  Response  for  25  Cases  with  Fourth  Degree  Actuators 
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Figure  G.65.  Bank  Angle  (<j>)  for  25  Cases  with  Fourth  Degree  Actuators 
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Figure  G.66.  Yaw  Rate  Response  for  25  Cases  with  Fourth  Degree  Actuators 
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Roll  Rate  Command  The  roll  rate  command  has  a  magnitude  of  90  de¬ 
grees/second  for  0.5  seconds. 


Figure  G.71.  Pitch  Angle  (0)  for  25  Cases  with  Fourth  Degree  Actuators 
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Figure  G.72.  Roll  Rate  Response  for  25  Cases  with  Fourth  Degree  Actuators 


Figure  G.73.  Bank  Angle  (<f>)  for  25  Cases  with  Fourth  Degree  Actuators 


Figure  G.74.  Yaw  Rate  Response  for  25  Cases  with  Fourth  Degree  Actuators 
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Figure  G.76.  Normal  Acceleration  for  25  Cases  with  Fourth  Degree  Actuators 


G-46 


Figure  G.77.  Lateral  Acceleration  for  25  Cases  with  Fourth  Degree  Actuators 
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